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Abstract 



The main goal of the study presented in this thesis was to perform in-situ in- 
vestigations on deformation structures in plastically deformed polycrystalline 
copper at low degrees of tensile deformation (< 5%). Copper is taken as a 
model system for cell forming pure fee metals. 

A novel S5mchrotron-radiation based technique High Angular Resolution 3DXRD 
has been developed at the 1-ID beam-line at the Advanced Photon Source. The 
technique extents the 3DXRD approach, to 3D reciprocal space mapping with 
a resolution of w 1 • 10~''A~^ and allows for in-situ mapping of reflections from 
deeply-embedded individual grains in polycrystalline samples diiring tensile 
deformation. 

We have shown that the resulting 3D reciprocal space maps from tensile de- 
formed copper comprise a pronounced structure, consisting of bright sharp 
peaks superimposed on a cloud of enhanced intensity. Based on the integrated 
intensity, the width of the peaks, and spatial scanning experiments it is con- 
cluded that the individual peaks arise from individual dislocation-free regions 
(the subgrains) in the dislocation structure. The cloud is attributed to the dislo- 
cation rich walls. 

Samples deformed to 2% tensile strain were investigated under load, focusing 
on grains that have the tensile direction close to the (100) direction. It was 
found that the individual subgrains, on average, are subjected to a reduction 
of the elastic strain with respect to the mean elastic strain of the grain. The 
walls are equivalently subjected to an increased elastic strain. The distribution 
of the elastic strains between the individual subgrains is found to be wider 
than the distribution of strains within the individual subgrains. The average 
properties are consistent with a composite t5rpe of model. The details, however, 
show that present imderstanding of asymmetrical line broadening have to be 
reconsidered. 

Based on continuous deformation experiments, it is found that the dislocation 
patterning takes place during the deformation, and that a subgrain structure 
appears from the moment where plastic deformation is detected. By investi- 
gating samples imder stress relaxation conditions, and unloading, it is found 
that the overall dislocation structure only depends on the maximum obtained 
flow stress. However, some changes in orientation and internal strain distribu- 
tion between the subgrains were observed after the unloading. 

An in-situ stepwise straining experiment of a pre-deformed sample was per- 
formed, allowing for investigation of individual subgrains during straining. 
The result indicates that the cell refinement process generally does not take 
place through simple subgrain breakups. Surprisingly, the dislocation struc- 
ture shows intermittent behavior, with subgrains appearing and disappearing 
with increasing strain, suggesting a dynamical development of the structure. 



ii 



Preface 



This thesis is submitted in partial fulfilment of the requirements for obtaining 
the Ph.D. degree in physics at Roskilde University (RUC). 

The research presented here was carried out within the Center for Fundamen- 
tal Research: Metal Structures in Four Dimensions (M4D), at Rise National 
Laboratory. The studies were conducted during the period from January 2004 
to December 2006. 

My supervisors have been Tage Emil Christensen (RUC), Hennrng Friis Poulsen 
(M4D) and Wolfgang Pantleon (M4D), who all are thanked for their help and 
support throughout the project. 

The experimental work was conducted at the 1-ID beam-lrne of the Advanced 
Photon Source (APS) at the Argonne National Laboratory, USA. The experi- 
ments would have been impossible without the close collaboration with our 
main contact at APS, Ulrich Lienert. He is warmly thanked for all his very 
valuable help, and for teaching me so much about experimental X-ray physics 
through all the long night shifts we shared. 

Beside U. Lienert I had the opportunity to meet and work together with a 
number of very pleasant people at the APS: John Aimer, Sarvjit D. Shastri, 
AU Mashayekhi, Joel Bernier, and Dean Haeffner. I am very grateful for all 
their help, without which we could never have performed the experiments 
presented, and for making all my visits to APS very pleasant despite the work- 
load. 

Carsten Gundlach (M4D), Henning Osholm Serensen (M4D), and Matthias 
Prinz (Freiberg University of Mining and Technology) are acknowledged for 
participating in beam times, and M. Prinz furthermore for doing some of the 
"data mining" for the Grain III dataset presented in section 14^21 

AU electron microscopy investigations conducted in connection with this study 
have been done by the "microscopy and sample preparation experts" in the 
M4D group: Qrngfeng Xing, Xiaoxu Huang, Gitte Christiansen, Preben Olesen, 
Hebner Nilson, and Guilin Wu. General preparation of the sample material 
was done by Palle Nielsen and Lars Lorentzen. They are all thanked for their 
very valuable help in producing and characterizing the sample material. 



iii 



iv 



Brian Ralph (School of Engineering and Design, Brunei University, UK) and 
Rasmus Brauner Godiksen (M4D) took time for reading and commenting this 
thesis, for which I am deeply grateful. 

The whole M4D group is thanked for making these 3 years very pleasant; espe- 
cially the "227 office"; Tine Knudsen, Kristoffer Haldrup and Rasmus Brauner 
Godiksen with whom I have sheared most of my time at Riso. 

Finally I would like to thank my long time colleague and friend Kristine Niss, 
and my wife Bodil Hjort Mynster for all the moral support that got me through 
the time as a Ph.D. student. 



Bo Jakobsen 
Ris0, December 2006 



Comment for arXiv/cond-mat version of the thesis 

The originally submitted thesis included the six papers quoted as Paper I - 
Paper VI (see page|99) as an appendix. Due to copyright issues the appendix 
has been omitted in this version, and minor changes has been applied to the 
main text accordingly. The thesis is meant to be self contained, but it is strongly 
advisable to also acquire the papers. 

The thesis was successfully defended Marts 2007. 



Bo Jakobsen 
Roskilde University, August 2007 



This work was supported by the Danish National Research Foundation and 
the Danish Natural Science Research Council. 



Use of the Advanced Photon Source was supported by the U. S. Department of 
Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. 
W-31-109-ENG-38. 



V 



A number of papers have been written in connection with this Ph.D. project, 
they are listed below. 

B. Jakobsen, H. F. Poulsen, U. Lienert, J. Bemier, C. Gundlach, and 
W. Pantleon. Stability of dislocation structiires in Cu towards strain 
relaxation. In preparation. 

B. Jakobsen, U. Lienert, J. Akner, H. F. Poulsen, and W. Pantleon. 

Direct observation of strain in bulk subgrains and dislocation walls 
by high angular resolution 3DXRD. Materials Science and Engineer- 
ing: A, 2007. Article in Press. doi:10.1016/j.msea.2006.12.168 

B. Jakobsen, H. F. Poulsen, U. Lienert, X. Huang, and W. Pantleon. 

Investigation of the deformation structure in an aluminium mag- 
nesium alloy by high angular resolution three-dimensional X-ray 
diffraction. Serif ta Materialia, 56:769-772, 2007. 

B. Jakobsen, U. Lienert, J. Aimer, W. Pantleon, and H. F. Poulsen. 
Properties and dynamics of bulk subgrains probed in-situ using a 
novel x-ray diffraction method. Materials Science Forum, 550:613- 
618, 2007. 

B. Jakobsen, H. F. Poulsen, U. Lienert, and W. Pantleon. Direct de- 
termination of elastic strains and dislocation densities in individual 
subgrains in deformation structures. Acta Materialia, 55:3421-3430, 
2007. 

U. Lienert, J. Aimer, B. Jakobsen, W. Pantleon, H. F. Poulsen, D. Hen- 
nessy, C. Xiao, and R. M. Sute. 3-dimensional characterization of 
polycrystalline bulk materials using high-energy synchrotron radi- 
ation. Materials Science Forum, 539-543:2353-2358,2007. 

W. Pantleon, B. Jakobsen, U. Lienert, J. Aimer, C. Gimdlach, and 
H. F. Poulsen. In-situ observation of individual subgrains by 3DXRD 

during deformation and recovery. In Proceedings of PLASTICITY 
'06: The Twelfth International Symposium on Plasticity and its Current 
Applications, pages 664-666, 2006. 

B. Jakobsen, H. F. Poulsen, U. Lienert, J. Akner, S. D. Shastri, H. O. 
Sorensen, C. Gundlach, and W. Pantleon. Formation and subdivi- 
sion of deformation structures. Science, 312:889-892, 2006. 

H. O. Serensen, B. Jakobsen, E. Knudsen, E. M. Lauridsen, S. F. 
Nielsen, H. F. Poulsen, S. Schmidt, G. Winther, and L. Margulies. 
Mapping grains and their d5mamics in three dimensions. Nuclear 
Instruments and Methods in Physics Research Section B, 246:232-237, 
2006. 

U. Lienert, J. Aimer, B. Jakobsen, H. F. Poulsen, and W. Pantleon. 
Observation of dislocation structure evolution by analysis of X-ray 
peak profiles from individual bulk grains. In Proceedings of the 25th 
Ris0 International Symposium on Materials Science: Evolution of Defor- 
mation Microstructures in 3D, pages 417-422, 2004. 



vi 



Dansk Resume 



Hovedformalet med det ph.d.-studie, som praesenteres i denne afhandling, er 
at foretage in-situ undersegelser af dislokationsstrukturer i plastisk deformeret 
kobber ved sma deformationsgrader (< 5%). Kobber skal i denne sammen- 
haeng ses som et modelmateriale for de rene fee metaller, hvor dislokationerne 
arrangerer sig i en celle-struktur. 

Vi har udviklet en synkrotronbaseret rontgenteknik (High Angular Resolu- 
tion 3DXRD), som tilf0jer hojoplost 3D kortlaegning af det reciprokke rum til 
raekken af 3D rontgendiffraktionsteknikker (3DXRD teknikker). Udviklingen 
af teknikken er foregaet pa 1-ID beam-linjen pa s5mkrotronen Advanced Pho- 
ton Source ved Argonne National Laboratory, USA. Metoden gor det muligt at 
unders0ge forbreddede Bragg refleksioner fra dybtliggende individuelle korn 
i en polykrystal, og det med en oplosning pa w 1 • 10~^A~^. Videre er tidso- 
plosningen god nok til, at man er i stand til at folge strukturudviklingen in-situ 
under deformation. 

Vi har fundet, at sadanne 3D kort over Bragg refleksioner fra plastisk deformeret 
kobber indeholder en udtalt struktur bestaende af skarpe toppe med hoj inten- 
sitet overlejret pa en sky af forholdsvis lav intensitet. Den integrerede intensitet 
i og bredden af disse toppe samt den rumlige fordeling af materialet, som giv- 
er anledning til disse, er blevet analyseret. Pa den baggrund konkluderer vi, at 
toppene er diffraktionssignalet fra de individuelle dislokationsfrie omrader i 
strukturen (under kornene). Tilsvarende konkluderer vi, at skyen stammer fra 
de dislokationsfyldte vaegge, som separerer under kornene. 

Underkornene i plastisk deformeret metal er traditionelt blevet undersogt med 
transmissionselektronmikroskopi, som er en destruktiv teknik, eller med k- 
lassiske rontgenteknikker, som midler over mange underkorn. "High Angu- 
lar Resolution 3DXRD" teknikken giver mulighed for direkte og med rimelig 
tidsopl0sning at undersoge egenskabeme af imderkom dybt inde i et korn. 

Den interne fordeling af elastisk tojning blev undersogt i prover deformeret til 
2% forlaengelse. Specielt fokuserede vi pa korn, hvor traekretningen er taet pa 
en (100) krystallografisk retning. Det blev fundet, at underkornene i gennem- 
snit er udsat for en reduceret elastisk tojning i forhold til den gennemsnitlige 
elastiske tojning i komet. Det blev tilsvarende fimdet, at vaeggene er udsat for 
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en forhojet elastisk tajning. Analysen viste desuden, at fordelingen af elastisk 
tejning mellem underkomene er bredere end fordelingen af elastisk tejning, 
som findes internt i de enkelte iinderkorn. 

Disse resultater peger pa, at komposit-modeller for dislokationsstrukturen beskriv- 
er den gennemsnitlige fordeling af elastisk tajning korrekt. Resultaterne vis- 
er imidlertid ogsa, at den nuvaerende forstaelse af asymmetrisk-forbreddede 
rontgen-linjeprofiler skal revideres, da man i eksisterende analysemetoder an- 
tager, at tojningsfordelingen mellem imderkornene er meget smallere end den 
interne tejningsfordeling. 

Et andet resultat af malingeme er, at taetheden af dislokationer i underkomene 
er meget lille. Teefheden af redimdante dislokationer er mindre end 12 • lO^^m"^. 
Vi har desuden fimdet, at teknikken kan vaere folsom for ned til en uparret dis- 
lokation. 

Pr0ver er ogsa blevet unders0gt under kontinuerlig deformation bade fra fuldt 
udglodet tilstand og fra pre-deformeret tilstand. Ved at folge en Bragg refleks 
fra et korn under deformation fandtes, at dislokationsstrukturen opstar og ud- 
vikler sig kontinuerligt under deformationen, og at underkomsstrukturerne 
eksisterer fra det ojeblik, hvor den plastiske deformation er detekterbar. Ved 
tilsvarende undersogelser, hvor spaendinger blev relakseret og proven aflas- 
tet, fandtes, at den overordnede dislokationsstruktur udelukkende afhaenger 
af den maximalt opnaede flydespaending. Under aflastning sas dog eendringer 
i tojnings- og orienteringsfordelingen mellem underkomene. 

En pre-deformeret prove blev undersogt under trinvis deformation. Sadanne 
forsog tillader, at man undersoger udviklingen af individuelle underkorn som 
funktion af tojning. Resultateme indikerer, at underkom ikke deler sig gen- 
nem en simpel opbygning af nye vaegge. De ser ud til at opsta og forsvinde 
d5mamisk, mens proven deformeres. Hvis dette resultat kan eftervises ved hoj- 
ere deformationsgrader, kan det f.eks. give en forklaring pa, hvordan disloka- 
tionsveegge opretholder en foretrukken orientering under deformation. 

De opnaede resultater kan forhabentlig inspirere til nye modeller for deforma- 
tionshaerdning og strukturformation i metaller. 
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List of commonly used notations 



q Scattering vector. Normally expressed in the defined q-space 

coordinate system 

6 Scattering angle 

60 Scattering angle for an undeformed perfect sample 

f] Azimuthal angle 

to, <p, X Available rotations on the setup 

Alo Rocking interval 

(}x, (jz, (}z Basic vectors for the reciprocal space coordinate system 

X, y, z Basic vectors for the laboratory coordinate system 

(fl, fj) Plane spanned by the two vectors a and h 

\uvw\ Lattice vector in crystallographic coordinates 

(uvw) Lattice vector family 

(hkl) Lattice plane with Miller index hkl 

{hkl} Lattice plane family with Miller index hkl 

hkl Reflection (or reflection family) from lattice plane (hkl) (or lat- 
tice plane family {hkl}) 

E Energy of the X-ray beam 

A Wavelength of the X-ray beam. 

L Horizontal sample-to-detector distance. 

ai, a2, 03 Crystal lattice basic vectors 

flp fl2' ^3 Reciprocal lattice basic vectors 




Sketch of the setup, defining axes in 
real and reciprocal space, rotation an- 
gles and scattering angles. For full fig- 
ure text see pagel4T] 
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Chapter 1 

Introduction 



This thesis deals with the fundamental properties of dislocation structures that 
evolve in metals as they are plastically deformed. 

There are, at least, two reasons why such structures are interesting. Firstly, they 
are an example of natural pattern formation which has a range of unexplained 
phenomena associated. Secondly the structures have an influence on the prop- 
erties of metals, and are hence interesting from an applied perspective. The 
first point is what is important to me, and I will not in this thesis discuss the 
consequences of the results in an applied framework. 

Two main issues have been investigated: 

• What are the static properties of the dislocation structures under load? 

• What are the dynamics of the structures during deformation? 

The main focus has been on pure fee metals, and the model material of choice is 
polycrystalline copper (Cu), which has been investigated in great details in the 
past. The deformation mode has been restricted to tensile deformation. Mainly 
small plastic deformations (in the range of — 5%) have been considered, as the 
focus is on the creation of deformation structures and their properties in the 
initial phase of structure formation. Focus has further been on grains where 
the tensile axis is close to the crystallographic (100) directioi^^- 

A major part of my work has been devoted to the development of a novel tech- 
nique: "High Angular resolution 3DXRD". The method is based on 3D reciprocal 
space mapping of individual reflections from individual bulk grains in poly- 
crystalline samples. By means of a setup developed at the 1-ID beam line of 
the Advanced Photon Source (APS) at the Argonne National Laboratory, USA, 

^Two reasons exist for this particular choice, firstly such grains are easy to locate with the used 
technique, and secondly it allows for comparison with some classical X-ray investigations. 
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CHAPTER 1. INTRODUCTION 



such 3D maps can be acquired in-situ and reasonably fast. These maps turn 
out to provide access to direct information on the deformation structure. The 
contents of the thesis reflects this experimental development. 

This thesis is divided into four major chapters: 

Chapter 1: Beside this short introduction, the present chapter includes an in- 
troduction to deformation structures in metals, the techniques used to 
investigating them, and the questions that I have touched upon in this 
study. This is followed by a brief overview of my Ph.D. project. 

Chapter 2: Gives an overview of the background of my work. This includes 
diffraction theory for deformed and undeformed metals, and classical X- 
ray methods for investigating such diffraction signals (section [2 . 1 1 to IZ4t . 
Section 12.51 discusses some of the recent developments in s5mchrotron- 
based techniques, two of which this work is based on. 

Chapter 3: Gives a thorough description of High Angular resolution 3DXRD. 
The chapter provides a detailed description of the setup, data analysis 
and arguments for interpretation. Finally the technique is compared with 
other techniques. 

Chapter 4: The major scientific re sults of my study is presented in a number 
of papers ( Paper ll4Paper v"l) . Chapter 4 gives an overview of the results. 



connecting results presented in different publications and presents some 
yet unpublished results. 

Conclusions and outlook are finally presented in chapter |5l 

References to the major publications are throughout this thesis designated as 



Paper IHPaper Vll, accordingly to the list presented first in the bibliography. 



1.1 Deformation structures in metals 

The plastic deformation of metals is carried by the movement of line defects, 
dislocations, through the crystal. 

During the propagation of the dislocations, some of them will be trapped in 
the crystal. The trapping can be due to e.g. foreign particles in the crystals, but 
stems mainly from interaction between individual dislocations. 

It is well known that the dislocations stored in a crystal have a tendency to self 
organize into ordered structures, known as dislocation structures or deformation 
structures. The morphology of the structures formed depends on the material 
investigated, the mode of deformation, and the degree of deformation. 
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Figure 1.1: Left) Atomic positions around an edge dislocation in a simple cu- 
bic (sc) crystal. The positions are calculated from the continuum displace- 
ment field for one edge disloca tion in an infinite crystal (as e.g. found in 
IWeertman and Weertmanl.ll964|l ). The position of the dislocation is indicated 
by the red symbol. A Burgers circuit around the dislocation is shown (blue 
line). Elastic properties of copper were used in the calculation, and the Burg- 
ers vector set to one inter-atomic distance. Right) Equivalent Burgers circuit 
in a perfect crystal. The corresponding Burgers vector is indicated by the red 
arrow. 



An introduction to dislocations and dislocation structures will be presented 
in the following. The focus will be on the material of choice, copper, but the 
results presented are general for a large cl ass of pure fee materials (for a general 
review of deformation structures see e.g. iHughes and Hansenl,l200l ). 



1.1.1 Dislocations 

A vast amount of literature exists on dislocation theory including a number of 
very good text-books such as IWeertman and Weertmanl,[l9641 and IHuU and Bacon , 
1984|| . for a general introduction to the concept please refer to one of these. 

The left part of figure lTD shows the distorted atomic lattice around one disloca- 
tion. A Burgers circuit is drawn around the dislocation. In the right part of the 
figure the corresponding circuit is shown in a perfect crystal, and the resulting 
closure failure shown; this is the Burgers vector of the dislocation. Dislocations 
are beside the Burgers vector characterized by their direction. Dislocations are 
separated into edge, screw and mixed dislocations. Edge and screw disloca- 
tions have Burgers vector perpendicular and parallel to the direction of the 
dislocation respectively, and a mixed dislocation an intermediate angle. It is in 
figure lLTl seen that the lattice is rather perfect far from the dislocation, and that 
a shear deformation is created if the dislocation is moved completely through 
the lattice. 
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D 



Figure 1.2: A) Symmetrical tilt boundary with dislocation spacing D and tilt 
angle /5. B) dislocation dipole. 

Dislocations give rise to displacement of the crystalline material around them. 
The description of a dislocation is normally divided into a part describing the 
"core" of the dislocation, that is the atoms very close to the dislocation, and a 
part describing the displacement of the crystal further away from the disloca- 
tion. 

Elasticity theory can be applied away from the dislocation core. Such elastic 
descriptions which gives the displacement, stress, and strain fields aroimd the 
dislocation, exists for multiple dislocation types and boundary conditions (see 
e.g. [Leibfried and Liicke, 1 9491 and the general discussion in e.g. 



An interesting and general feature of the stress / strain fields of a dislocation is 
that they are long-ranged. They go to zero as 1 / r (r being the distance from the 
dislocation). This has the consequence that the elastic energy of a dislocation 
generally increases logarithmically with the size of the crystal, and therefore 
diverges. 

Dislocations can be arranged in a number of stable configurations, with lower 
energy than a single dislocation. Figure FOl show two examples of such struc- 
tures, the symmetrical tilt boundary and the dislocation dipole. The stress field of 
the symmetrical tilt boundary goes down as fast as e~' , and the dipole field 
as l/r'^. The tilt boundary has the property that the crystal parts on each side 
are rotated with respect to each other. The angular difference between the two 
sides, fi is given as j6 = / D with b the length of the Burgers vector and D the 
spacing between the dislocations. 

An interesting configuration of dislocations is a spatially random distribution 
of dislocations, having an equal number of dislocations of opposite Burgers 
vector. It can be shown that the elastic energy for such a distribution a lso di- 
verges logarithmically with the size of the crystal lWilkensl,ll970bl,ll984ll . 

Hence in order to lower the energy the dislocations stored in a crystal have to 
be arranged in ordered structures. 
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1.1.2 Phenomenology of dislocation structures in copper 



The deformation structures in copper after tensile deformation have been in- 
vestigated in great detail by electron microscopy over the last 40 years or t 
both singl e crystals (to name a few: [Essmann , 'l963; St eedsl Il966t iGottler , 
[l973tiKawasaki and TakeuchiLll980Hwilkens et aZ.l ll987]) and polycrystals (e.g. 
I Essmann et flZ.LIl968l;lHuangLll998ll V 



It is generally found that the dislocations self organize into regions of rela- 
tively low dislocation density, knows as cell interiors or subgrains, separated 
by regions of much higher dislocation density, known as cell walls or subgrain 
boundaries. An example, in the form of an electron micrograph, of such a struc- 
ture is shown in figure 11.31 The cell walls seen in the micrograph are rather 
loose, as the sample was only strained to 2% tensile deformation. At higher 
plastic deformation the walls becomes sharper. 

This separation into cells and cell walls are general for a large class of ma- 
terials known as cell-forming metals (e.g. Al an d Ni). For a general review 
of deformation structures see e.g. LHansen and Jensen , 19991: iHansenl, I2OOII: 
Hughes and HansenLlzOOl . 



The morphology of the dislocation structure does, as mentioned earlier, de- 
pend on the deformation mode, but even with the same deformation mode 
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large differences may exist. In the case of unidirectional tensile deformation it 
is found that the morphology depends on the direction of the tensile axis with 
respect to the crystallographic orientation of the crystal, with equivalent be 



havior for sin gle cr ystals and grains in polycrystals. See [Huans an d Hansen , 
19971: Ih uangi, Il998|l for work on polycrystals and references therein for work 



on single crystals. 

In the general case the structure is hierarchal, consisting of what is known as 
cell blocks which again is separated into cells. The cell blocks are separated by 
boundaries of a rather high misorientatio n, whereas the cells are separated by 
boundaries of low misorientation [Kuhlma nn-Wilsdorf and Han sen, 1991]. 

In my work, I have focused on crystals which have a (100) direction close to the 
tensile axis. The morphology in this case consists of eq uiaxed cells se parated 
by cell boimdarie^as the one seen in figure [Ol (see e.g. iG6ttleij,[l973ll '). 

A very pronounced feature of the dislocation structure is that the length scale 
decreases with increasing plastic strain, a phen omenon known as cell refine- 
ment (see e.g. the classic work by lGottleil lll973ll ). The misorient ation between 



the cells furthermore increases with increasing plastic strain (e.g. HHughes et al 
19971: |Pantleo^.l2002ll V 



The underlying principles controlling the structural formation have been a 
matter of debate over the last many years. Two main ideas exist on why struc- 
tures form: lowering of energy (Low Energy Dislocation Structure (LEDS)), 
and self organization. In the case of LEDS theory it is argued that the dislo- 
cations will form structures that among the available confo rmations lead to 



a minimization of the free energy (e.g. IIKuhlmann-Wilsdorii 1200111 ). Self or- 
ganization theories are based on the general observation that driven systems 
far fr om equilibrium have a natural tendency to form structures (e.g. [SeegeJ 



1988(1 ). F or a detailed discussion of different models see the comprehensive 



review by lKubinllll992ll . 

However, it seems to be generally accepted t hat dis locations n eed to be mo- 
bile in three dimensions [ Kuhlmann-Wilsdorf et al 1 IT994: Madec et al\. 120021: 



Hughes and HansenI, |2004| for structural formation to exist. This has the con- 



sequence that properties such as stacking fault energy and alloying elements 
change the structure formation ability of a metal. 



1.1.3 Techniques for investigating dislocation structures 

Classically two main methods are used for the investigation of dislocation 
structures; transmission electron microscopy (TEM) and X-ray line profile anal- 
ysis. 



hierarchal structure als o exists for t his crystal orientation, consisting of the cells and what 
is known as supercells tWilkens et g/.L[l98^ . 
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Transmission electron microscopy 



Transmission electron microscopy is by far the most commonly used technique 
for investigating dislocation structures. The technique gives very informative 
real-space images of the dislocation structures but has some disadvantages, 
mainly related to the fact that thin films have to be prepared for investigation: 



Care has to be taken to hinder relaxation of dislocation structures and in- 
ternal strain distribu tion during sample preparation . 

See e.g. [ E ssma n^ 1963 : Young and Sherrill , 19671: Crump and Young , 
19681: lMughrabil,ll97lF 



• It is impo ssible to investigate bu lk samples under tensile deformation. 
See e.g. iMyshlyaev et ali 1978 | ] for an example of creep investigations 
performed on thin films, and iMartin a nd Kubin, 1978] for a discussion 
of the limitation of in-situ studies by high voltage TEM. 



Traditional X-ray techniques 



Alternately; characterization can be performed with X-ray techniques , espe 
dally line profile analysis and rocking curve analysis (see e.g . jwilkens, 1970, 
Wilkens et flllll987tlm"rren and Averbachl |l950':'Unga r et fl/T[T984aciKrivogla: 



199611 '). Further detains on some of these techniques are provided in section 



23A\ 



The methods will give information on average properties such as strain distri- 
bution, dislocation density and orientation distribution. The main advantage is 
that the techniques are non-destructive. The techniques are on the other hand 
highly indirect as: 



• The results represent the average over many structural elements (sub- 
grains and in the case of polycrystalline samples also grains). These el- 
ements will have one common direction of the lattice plane normal, but 
have different crystallographic orientations in the sample and different 
neighboring environments. 

• Models are needed to interpret the results. The use of such models in- 
troduces additional assumptions which often can be hard to verify by 
independent experiments. 



Recent X-ray methods 



The availability of s5mchrotron radiation has lead to the development of new 
methods. Most interesting for the present work are the 3D X-ray diffraction 
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3DXRD) mi croscope f Poulsenl l2004ll , and the 3D X-ray crystal microscope 



Larson et a?.Ll2002i1 (described in section |2311 and |2.5.2t . 
Both techniques do however have some limitations: 

• The 3DXRD microscope does not have the spatial resolution for direct 
observation of the deformation microstructure. 

• The 3D X-ray crystal microscope is based on a spatial scanning technique, 
which limits the volume that can be mapped in a feasible time. 



1.1.4 Open questions 

Fundamental questions on deformation structures exist even though they have 
been investigated intensively over the last many years. During my work I have 
touched upon a few of such questions: 

• What is the elastic strain in the subgrains and subgrain walls? 
(section 112) 

• What is the dislocation density in the subgrains? 
(section |43) 

• When do the dislocation structures form? 
(section HXJl 

• What is the stability of dislocation structures? 
(section 11331 

• How does the cell refinement take place? 
(section |43) 



1.2 Development of the project 

A close interlink exists between the questions dealt with, and the development 
of High Angular Resolution 3DXRD. 

The original plan for the Ph.D. project was to use "3DXRD peak shape anal- 
yse" (see section |2. 5.31 1 for in-situ investigations of the dynamical properties of 
dislocation structures. Such a study would have been rather simple from an 
experimental perspective, as the technique already existed. As for many other 
X-ray techniques it would have required a heavy use of models for interpreting 
the data. 

However, during one of the very first beam-time experiments we discovered 
that it might be possible to obtain diffraction signals from individual subgrains 
in a deformed structure. 
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Based on these first indications, with the experiments conducted on thin films 
and on polycrystalline aluminum samples, it was decided to use more beam- 
time on exploring these possibilities. In close collaboration with U. Lienert (oui 
local collaborator at APS) it was decided to change the X-ray optics completely 
with respect to what had been used before. This new X-ray optics setup al- 
lowed for reasonable data acquisition tunes and high angular resolution at the 
same time. 

Based on the data acquired it was confirmed that we were able to obtain in-situ 
diffraction data directly from individual subgrains embedded in a bulk grain 
in a polycrystalline sample. 

With the technique established a range of questions came to mind which might 

be possible to investigate in a more direct way than previously had been pos- 
sible. This included the questions on subgrain d5mamics which was the initial 
goal, but also questions regarding the internal strain distribution in the grain 
and the consequences for line profile analysis. 
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Chapter 2 



Background 



The relevant basic diffraction theory and general experimental methods will be 
briefly reviewed in this chapter. 

It should not be seen as a general introduction to diff raction, for thi s I refer to 
the la rge n umber of text books on this subject, such as iGuinieiill963ll , IWarrenl 
1969ll and llAls-Nielsen and McMorrowllZOOlll . 



2.1 Basic scattering theory 

General scattering theory will be described in the following section. The theory 
will be restricted to kinematic scattering in the elastic limit of monochromatic 
X-rays, as these are the relevant conditions for the present study. 

An object consisting of a number of atoms is shown in figure lZTl The scattering 
ability of the individual atoms is described by the atomic scattering factor, fj, 
and the position by the vector, rj. 

The object is illuminated by a plane wave monochromatic X-ray beam de- 
scribed by the wave vector ko. The scattered wave is observed at the point, 
O. This observation point is assumed to be far away from the object relative 
to the size of the object, therefore the scattered waves from the different atoms 
can be described by the same wave vector k. The length of the wave vector is 
preserved due to the assumption of elastic scattering, that is: 

\k\ = \ko\ = ^ (2.1) 
where A is the wavelength of the X-ray beam. 

The scattering property of such an object is described by the complex scatter- 
ing amplitude. A, which describes both the amplitude and the phase of the 
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Figure 2.1: A general ensemble of atoms. The scattering ability of the j'th atom 
is described by the atomic scattering factor for that atom, fj, and its position by 
the vector ry. ko is the wave vector of the incoming monochromatic beam, and 
k is the wave vector of the scattered beam as observed from the point O. It is 
assumed that the distance from the ensemble of atoms to the observer is much 
larger than the dimensions of the collection that is: |0 — r, | ^ k; ~ i'k\'^hj>k. 

observed scattered wave relative to the incoming wave. The phase difference, 
in the scattered waves, due to the different positions of the atoms can be found 
as: 

{k-ko)-rj = q-Tj, (2.2) 

with the scattering vector, q, defined as0 <J = fc — fco. 

The scattering amplitude from a collection of atoms can be written as: 

A{q)=Y.fii'iy'"'' (2-3) 

where fj{q) is the (/-dependent atomic scattering factor for atom 

However, X-ray detectors do not record both the phase and the amplitude of 
the scattered beam, but only the intensity, l{q), which is given as: 

l{q) = AA* = \A[q)\h (2.4) 
hence the phase information is lost. 



^The scattering vector is defined by some authors (e.g. lAls-Nielsen and McMorrowl llOOlll ) 
with the opposite sign. 
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2.2 Diffraction from a perfect crystal 

The crystal lattice 

The position of the atoms in a crystalline material is normally described by a 
lattice and a basis. 

A crystal lattice is characterized by the fact that it obeys certain translation 
s5Tnmetries. A 3D lattice can be described by three crystal lattice basis vec- 
tors, ai, fl2 and a^, which have the property that the lattice looks the same if 
translated by an integer number of any of these. 

The lattice is more formally described by vectors in the form: 

R„ = «iai +n2fl2 + «3«3/ (2.5) 
with n = («!, «2/ W3) all being integers. 

These vectors give the position of the unit cells of the crystal, the lattice points, 
each unit cell is populated by the same arrangement of atoms described by 
what is known as the basis. In total: 

lattice + basis = crystal structure. (2.6) 

The choice of crystal lattice basis vectors is not unique, nor is the basis. 

The basis can be described by vectors, ry, relative to the lattice points. The 
position of any given atom in a crystal can be given as: 

Rn^j = R„+ rj = nifli + «2a2 + «3«3 + rj, (2-7) 

for some n, j. 



Scattering amplitude 

The general formula for the scattering amplitude (equation |23) can in the case 
of a crystal be separated into two parts as: 

unit cell sum lattice sum 

Ml) = E/;(^)'^'' =E^(^)^E5^' (2-8) 

n,i 7 n 

where the "unit cell sum" is the sum over the atom configuration in the basis, 
and the "lattice sum" is over all lattice points. 
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The reciprocal space and lattice 

In the description of diffraction it turns out to be very useful to construct what 
is known as reciprocal space. 

The reciprocal space is spanned by the reciprocal basis vectors, a^, and fl|. 
These basis vectors are related to the crystal lattice basis vectors by: 

* 27r ^ 2n In 

flj = — ai X a3, — — fl3 x a\, aj, — — a\ x a-i, (2.") 

Vc Vc Vq 

with Vc — ai- {a2 X a^) the volume of the unit cell. It can be seen that the 
dimension of the reciprocal lattice vectors are reciprocal in length, hence the 
name. 

The two sets of basis vectors have the property that: 

a; • aj = 27tSij, (2.10) 
where Sij is the Kronecker delta. 

In the cubic case we have || aj, and | | = It is important to notice that 

the reciprocal space is tightly bound to the crystal. If the crystal is rotated so is 
the reciprocal space. 

The reciprocal basis vectors span, in a natural way, a lattice in reciprocal space, 
with a reciprocal lattice vector, G, given as 

Ghki = hal + kal + /a|, (2.11) 

with h, k, I integers. 

Reciprocal lattice vectors have the following properties, relating them to the 
imderlying crystal structure: 

• G;,j-/ is perpendicular to the lattice plane with Miller indices hkl. 

• \^hkl\ — where is the lattice spacing of the lattice planes with 
Miller indices hkl. 

The Laue condition 

The scattering vector can be described in coordinates of the reciprocal space in 
a natural way: 

q = qial + qja^ + qsa^, (2.12) 
with qi, q2 and q^ real dimensionless numbers. 
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The product • R„ in equation 12 . 8 1 then becomes: 

q -Rn =2n{niqi+n2q2 + n3q3), (2.13) 

using equation l2.10l 

In the case where <j is a reciprocal lattice vector, this sum reduces to an integer 
times 2n. The lattice sum in equation 12.81 then equals the number of lattice 
points (a large number) whereas it is of the order of unity in all other cases. 

This is the Laue condition for observation of X-ray diffraction^. 

q = G (2.14) 

When investigating the diffracted intensity as function of the scattering vector, 
a detectable signal is only obtained when the Laue condition is fulfilled. Such 
maxima are normally termed Bragg peaks or reflections. The measured inten- 
sity, in the individual reflection from a crystal, is determined by the unit cell 
sum for the crystal. 



2.3 Diffraction from real crystals 

The description in the previous section is based on a perfect infinite crystal, 
somewhat different from the crystals investigated in reality. The lattice sum 
will be non zero in some region around the theoretical reciprocal lattice point 
if the crystal has e.g. a finite size, a distribution of lattice spacings or a distri- 
bution of lattice plane orientations; the reflection is said to be broadened. 

The idea of many diffraction-based methods is to investigate such broadened 
reflections, and thereby obtain information on the material. 

As a crystal is deformed, the lattice structure becomes distorted. This can 
mainly happen in two ways: the lattice plane spacing can change, and the ori- 
entation of the lattice planes can change. A crystal will after plastic deforma- 
tion generally contain a distribution of lattice plane spacings and orientations, 
the latter sometimes termed the mosaic spread. 

From the two main properties of the reciprocal lattice vectors (as stated on pre- 
vious page) it is known that the length of a reciprocal lattice vector is related to 
the lattice plane spacing of the crystal, and that the orientation of the reciprocal 
lattice vector is related to the orientation of the crystal lattice planes. 

^The Laue condition can be shown to be equivalent to the Bragg condition for diffraction: 
2dsin(f?) = A where d is the lattice spacing for the relevant reflection and 8 the scattering angle. 
The advantage of the Laue formulation of the diffraction conditions in reciprocal space is that aU 
results and interpretations only depend on the underlying crystal structure. 
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A uniform straining of a crystal will hence lead to a radial shift, that is along 
the reciprocal lattice vector, of the reflectior0. An uniform rotation will equiv- 
alently lead to an azimuthal shift, that is perpendicular to the reciprocal lattice 
vector, of the reflection. 

A distribution of lattice plane spacings (equivalently elastic strains) will give 
rise to a broadening of the intensity distribution in the radial direction. The dis- 
tribution of orientations of lattice planes, will on the other han d give a broad 



ening of the intensity distribution in the azimuthal directions. llWilkenslll984 1 



The broadening is easy to understand if the crystal is thought of as consisting 
of a number of incoherently scattering domains, each with some strain and ori- 
entation. In this case one can think of the broadened reflection as being the 
simple superposition, in intensity, of a large number of reflections. The peak 
shape will in such a case be related, in a simple way, to the distribution of 
strains and orientations among the domains. However, a plastically deformed 
metal can in the general case not be divided into such incoherently scatter- 
ing domains with a well-defined strain and orientation, hence more elaborate 
models are needed taking into account the full distributions. 



2.3.1 Quantitative analysis of broadened reflections 

A few quantitative results regarding the broadening of reflections will be dis- 
cussed in the following. 



Broadening due to the finite size of the crystals 

Beside broadening due to strain and orientation the reflections will also be 
broadened because of the finite site of the crystals. 

Follow ing the derivation by lKrivoglazI |l996|l (originally formulated by lEwaldl 
| l940|| ) the general expression for the scattering amplitude (equation 12 8t can 
be reformulated as: 

Mcj) = EM^JV'"' I Y'"{r>{r)e"^-'dr, (2.15) 
/■ 

where y°° describes the position of unit cells in a infinite large crystal, and s(r) 
is a function which is 1 inside the actual crystal and outside. 

It is shown that the shape of the intensity distribution close to a reciprocal 
lattice point, G, is given by: 

W('7)'^ |s(«7-G)|2, (2.16) 

where s is the Fourier transform of s, and the unit cell sum has been neglected. 

^The change in length of the ij-vector from a reference length qo is for small changes, A(j, 
directly related to the strain in the crystal hd/ as hd/ dg = — Aq / qg to first order. 
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In the case of an infinite crystal we have s = 1 for all space and s{x) = S{x), 
(with S{x) being the Dirac delta function), equation 12.161 hence reduces to the 
exact Laue condition as expected. 

The width of the size-broadened peak, Aq, in some direction in reciprocal space, 
is related to the real-space length scale, /, of the crystal in the same real-space 
direction by: 

Acj = k^. (2.17) 

where k is the Scherrer constant, which is related to the precise measure of the 
peak width, and t he shape of the crystal (for a review of Scherrer constant in 
different cases see llLangf ord and Wilson , 1978 1). The Scherrer is generally not 



far from unity, it may for example be shown that k = 0.88 for a box shaped 
crystal and a width-measure of full width at half maximum. 



Strain broadening 

It is general for classic formulations of strain broadening that they relate the 
azimuthally integrated radial peak profile of the reflection to some description of 
the strain distribution in the crystal. The reason for studying such integrated 
radial peak profiles is that they can easily be obtained from single crystals and 
powders with conventional diffracto meters. 

The derivations below follow Warren and Averbach as described in 
I Warren and Averbachl il950i1 and [Warren. . 1969.] . 



The deformation is assumed to be so smooth that the individual unit cells in 
the deformed crystal are equal (hence it makes sense to talk about a unit cell 
sum). The unit cell sum will in the following, without loss of generality be set 
tol. 

The general equation for the diffracted intensity (equation l2.4t can be rewritten 
as: 

j(^) =^ ^^^e^-y-lKn-R,,,). (2.18) 



The position of the imit cells in a distorted crystal can be described by: 

R„ = + n2fl2 + W3«3 + (^l(n)fli + ^2(w)a2 + «^3(«)«3)/ (2.19) 

where S{n) = {Si{n),S2{n),S5{n)) is a small perturbation of the position of the 
unit cells. 

The result will be limited to the case of the radial intensity profile near an OO/q 
reflection, and the tj-vector is therefore written as: 

(7 = (jlflj +r?2fl2 + Co + t?3)«3/ (2.20) 
with (ji, q2 and small quantities. 
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Equation 12.181 can to a first order in the small quantities {S and (?i,^2/'?3) be 
written as: 

J = ^^g27r/({«i-n;)<?i + (n2-«9<72+(n3-«3)0+'o('53(«)-'53(«'))^ (2.21) 
« n' 

with I = Iq + qo,. 

By integrating this equation in qi and q2 over the full peak one obtains the 
integrated radial line profile: 

J = ^^^^e27ri{(n3-n^);)g27r/[{^3(«3)-^3(«^))/o]„j,„2 ^ (2.22) 
"1 "2 "3 «^ 

For a fixed ^3, the sums can be interpreted as a sum over all pairs of cells 
which are in the same column (along the aj, direction) and have a distance 
of («3 — n3)|fl3|. Now let N„ be the number of such pairs, and introduce the 
abbreviations n = — n'^ and Sn = ^^{n^) — S^{n'^). This reduces equation 
Eiato: 

^= E N„ (e^^'MA e27r//n (2.23) 



with () being the mean value over the entire crystal. 

From this rather long exercise, it can be seen that the radial intensity profile in 
a natural way can be written as a Fourier sum. The N„ term is related to the 
size broadening and will be ignored, hence the n'th term in the Fourier sum of 
the strain broadened profiles, A^^^^^{n), is given as: 

^strainjij^-j ^ ^^IniloS,,^ = ^^Imlone,,^^ (2.24) 

where Cn = S„ / n is the strain taken over a distance of « | fl3 1 in the direction of 
as- 

By approximations of this equation and models for the strain distribution it 
is possible to obtain analytical relations between the differential strain or e.g. 
dislocation dist ributions and the peak profile. Examples of such are th e clas- 
sical studies by lWarren and Averbach [,1950.1 , , Wilkens. l,1970a.d and .Krivog laz 
Jl99i. 

A simple relation exists bet ween the integral width , /3 of the intensity profile 
and the Fourier coefficients I van Berkum et al\ 1999 1: 

/oo 
A'''^'''{n)dn, (2.25) 
-00 

where the Fourier coefficients have been generalized to a continuous variable. 
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If the differential strain distribution, eg, in the material is assumed to be Gaus- 
sian it can be shown that: 

/5 = ^V^<e§>/o, (2.26) 

it is, as expected, seen that a close relation exists between the peak width 
relative to the peak position (/S/Zq) and the width of the strain distribution 



Asymmetric line broadening 

Beside simple broadening of the integrated radial peak profiles, it has been 
found that the profiles from plastically deformed metals in some cases show a 
pronounced asymmetry. This was first observed by Ungar during the study of 
tensile deformed single crystals of copper [ ,Ungar et a/.„,1984ajl . 

The obse rved asymmetry has been rationalized o n the basis of the comp osite 



ine obse rved asymmetry nas peen rationalized o n tne basis ot tne comp ( 
model by lMughrabil lll983ll in lUngar et fl/.lll984allbl:lMughrabi et anil986ll . 



The basic idea of the composite model is that a deformation structure is re- 
garded as consisting of two parts; the interior of the cells, which are relatively 
soft and the walls which are relatively hard (due to the large dislocation den- 
sity). The left part of figure 12721 illustrates such a system. 

Internal stresses will exist, during and after, plastic deformation (as an example 
here taking simple tensile straining) due to the difference in yield strength in 
the two parts. A backwards stress (with respect to the applied external tensile 
stress) will exist in the interior of the cells, leading to a reduction in the total 
stress in the cells. A forward stress will in a similar manner exist in the walls. 

The system mainly studied by Ungar et al. is single crystals with the tensile 
axis in the [100] direction. This has two consequences; firstly it ensures that 
the dislocation structure is cell like (see section ri.l.2|l . Secondly it allows for 
easy investigation of the diffraction signal from lattice planes with lattice plane 
normal parallel to the tensile axis (known as the axial case) and with lattice 
plane normal perpendicular to the tensile axis (known as side cases). 

The lattice spacing observed in the axial case will be larger for the walls than 
for the interior due to the internal stress differences. The part of the radial peak 
profile arising from the cell interiors will hence be shifted to a higher radial q 
position, and the part from the walls to a lower. In the side case this will be 
reversed due to the cross contraction of the crystal. The peak from the walls 
will, at the same time, be very broad as a wide strain distribution exists here. 

It was shown I Ungar et al. , 1984all that the asymmetric peak profiles can be 



divided into two "well behaved" symmetrical parts, which can be interpreted 
as the signals from the two parts of the composite. The right part of figure 
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Figure 2.2: Left) Illustration of the composite model of a deformed cell form- 
ing metal under a tensile load, a is the applied stress , Aa,i, and Aac are the 
forwards and backwards stresses in the walls and cell interiors respectively. 
'?axial and <7side indicate the direction of the scattering vector for the axial and 
radial cases respectively. Right) Illustration of an asymmetric peak in the axial 
case (red line), and the decomposition into a part from the cell interior (blue 
line) and a part from the walls (green line). 

I2.2l illustrates this decomposition in the axial case. The two symmetrical peaks 
have then been treated by classical line broadening theory. 

This shows that by analysis of broadened reflections from plastically deformed 
metals it is possible, at least in the context of the composite model, to obtain 
detailed information on the deformation microstructure. The composite model 
will be discussed in relation to the present study in section l4!2l 



2.4 Experimental investigation of diffraction 

Three basic methods exist for investigating the diffraction signal from crys- 
talline samples: The rotation method, the Laue method, and the powder method. 

The first two of these methods has classically been used on single crystals, and 
the latter on powder samples. However, with the use of modern synchrotron 
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Figure 2.3: Typical setup for the rotation method, using a transmission geom- 
etry and an area detector. Data is obtained by rotation the sample over some 
angle around the y-axis while exposing. The scattering angle, 29, azimuthal 
angle, r], and rotation angle (p are defined. 

radiation it is possible to investigate single grains in polycrystalline samples 
(see section [2.5.11 and |2.5.2| |. The Laue method uses a polychromatic X-ray 
beam, whereas a monochromatic beam is used for the other two methods. 

The monochromatic techniques will be discussed in the following, as they are 
the basis for the developed technique. Reciprocal space mapping will further- 
more be discussed. 



2.4.1 The rotation method 

A typical setup for rnvestigatrng the diffraction signal by use of the rotation 
method in the transmission mode is shown in figure 1231 The sample is illumi- 
nated by a monochromatic beam and the diffracted beam is recorded on an area 
detector. Lattice planes that happen to fulfill the diffraction condition will give 
rise to a diffraction spot on the detector, and by rotating the sample, different 
reflections can be brought into the scattering condition. 

In the rotation method, data is obtained by rotating the sample with constant 
angular velocity over some angular range, A(p, around a fixed axis (on the fig- 
ure the i/-axis), while exposing. There are two reasons for this. Firstly, the re- 
flections are very close to delta shaped if the sample is a perfect crystal, hence it 
is very hard to align the sample precisely at the scattering condition. Secondly, 
it will lead to a integration over the full intensity in a reflection if it has some 
width in reciprocal space. 

Such a single exposure integrates over some part of reciprocal space. To sample 
a larger part of reciprocal space many exposures are taken at adjacent angular 
intervals. 
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The total inte grated diffracted intensity for one reflection E;,j-/, is given by 
I Warrenl[T969ll : 



^-y-' - ^ j vl sin(20;,«)|sin,| ^^'^'^ 

where co is the angular velocity used for the measurement, V£.j.ystai is the vol- 
ume of the scattering crystal, Zg is the input intensity (energy per time per area), 
P{^hkl' is the polarization factor which depends on the polarization of the X- 
ray beam, t] is the azimuthal angle for the reflection, 0/,,^/ is the scattering angle 
for the given reflection family, F;,j;; is the structure factor for the given reflec- 
tion family, Va is the volume of the unit cell and e, and c are fundamental 
constants: charge of the electron, mass of the electron and speed of light, re- 
spectively. The last term in the equation is normally called the Lorentz factor, 
and depends on the details of the setup. 

What is most interesting for the present study is that the intensity is linear in 
the diffracting volume. 



2.4.2 The powder method 

A powder sample consists of a large number of small (with respect to the beam 
size) crystallites. 

The crystallites have some distribution of orientations (known as the textur^oi 
the powder). Generally a large number of crystallites will fulfill the scattering 
condition when the sample is illuminated at any orientation. Data is normally 
obtained with a stationary sample, using a monochromatic beam. The signal 
on the detector consists of what is known as Debye-Scherrer rings. 

The integrated diffracted intensity per angular unit in one Debye-Scherrer ring 
from a powder with random orientation of the crystallites is given as: 

F -T f \ ^%o^derm\FMlf 1 At 

Epowder - ^0 j ^2 ^^^>^>^^>I^iMMi) 360^ ^ ^ ^ 

where Vp^^^dgj. is the volume of the powder that is illuminated, AHs exposure 
time, m is multiplicity of the reflection family, and other symbols have the same 
meaning as in equation l2.27l Equation |Z28] is derived under the assumption of 
an unpolarized beam (as from a conventional X-ray source), which simplifies 
the calculation. However, in the case of S5mchrotron radiation the beam is al- 
most linear polarized in the horizontal plane, which means that equation l2.28l 
in the general case has to be modified. As the general result is not needed for 
the present studies further discussion of this will be postponed to section |3.6.5l 



^Texture is a general term in material science, and is e.g. used about the distribution of orien- 
tations of grains in a polycrystalline sample. 
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Sample 





detector 



Figure 2.4: General schematic of a classic point detector-based diffraction 
setup, fco arid ^ are the incoming and outgoing wave vectors respectively, 6 is 
the scattering angle and Ao; the misalignment in the scattering plane between 
the ideal lattice plane normal and the scattering vector. 

2.4.3 Investigation of intensity distributions 

In section [23] it was discussed what can be learned from the shape of the in- 
dividual diffraction peaks, the reflections. The reflections will in general be a 
three dimensional, intensity distribution in reciprocal space, close to the corre- 
sponding theoretical reciprocal lattice point. 

3D intensity distributions have normally been investigated in a number of 
ways, such as: line profile analysis, rocking curve analysis, 2D reciprocal space 
mapping, and 3D reciprocal space mapping. Beside the 3D mapping, these 
techniques represent a projection onto either a line or plane of the full intensity 
distribution. 

The schematic layout of the classic experimental setups for these types of mea- 
surements are generally much alikqj Figure 12.41 shown a schematic of such a 
setup based on a point detector. In the following some of the different applica- 
tions will be discussed. 



^However many different variations exists, using e.g. a line detector instead of a point detector 
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Line profiles and rocking curves 

Classic integrated radial line profiles are obtained by integrating over the az- 
imuthal directions, and investigating the intensity as a function of the length 
of the scattering vector (keeping the angle between the sample surface and 
scattering vector constant). Such curves are normally obtained by 9,29 scans, 
where the sample and detector are rotated in steps of A9 and 2A9 respectively. 
It can be seen from figure |2]4| that such a scan will keep the direction of the scat- 
tering vector, q, constant with respect to the sample surface while changing the 
length. The integration over the azimuthal direction out of the scattering plane 
is normally obtained by having a rather large beam divergence in this direction. 
The in-plan azimuthal broadening is normally integrated over by rotating the 
sample over an angle, Alo, for each 9, 29 point. 

The intensity distribution in the azimuthal direction can be investigated by 
what are known as rocking curves. The sample is here rotated with a constant 
speed in small steps around an axis perpendicular to the scattering vector (as 
e.g. represented by the angle Aco on figure 12.411 while recording the intensity. 
An integration over the radial and the other azimuthal direction can be ob- 
tained by suitable combinations of beam divergence and energy spread in the 
X-ray beam used. 



Reciprocal space maps 

The two above techniques can easily be generalized to 2D and 3D dimension al 
reciprocal space maps (for a review of such techniques see e.g. | Fewster , 1997|l ^. 



Most common is a projection onto the plane of the radial direction and the in- 
plane azimuthal direction. Such a map can be obtained by gathering radial line 
profiles for different angles between the sample and the scattering vector. This 
corresponds to 9,29 scans with different Ao; in figure IZ4l This requires a beam 
with a narrow energy spread, and a low divergence in the scatteri ng plane 



(an e xample of a diffractometer enabling such maps can be found in iFewster , 

By limiting the divergence in both the in plane and out of plan directions, it is 
possible to obtain full 3D reci procal space maps by introd ucing a second rota- 
tion axis for the sample (e.g. iFewster and Andrew , 1995ll ). The problem with 



such techniques are that they are point-by-point in a 3D space, hence acquisi- 
tion time rises quickly with the resolution obtained. Furthermore the resolu- 
tion tends to be different in the three directions. 



2.5 Recent developments 

The Laue and rotation methods have normally been used on single crystals. 
However, with the availability of synchrotron radiation and X-ray optics they 
have been generalized to multi-grain samples. 
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These generalizations are based on the fact that the individual grains in a multi- 
grain sample can be investigated using the classic methods if the beam size, 
sample properties (such as thickness and grain size) and other experimental 
parameters are matched. 

In the following three methods will be discussed as they have special relevance 
for the problems investigated here. 

It should also be mentioned that a large number of experiments which nor- 
mally were performed on home sources can now be perfo rmed in improved 
versions using synchrotron radiation (a few examples are | Chang et ah . 1199.4 



iBiermann et d.Ll997: Murph y et al.lQOh Schafler et al.i\lQQ^ ). These improve 



ments are normally better time and spatial resolution due to the much larger 
flux. 



2.5.1 3DXRD microscopy 

The 3D X-ray diffraction (3DXRD) method is a monochromatic S5mchrotron- 
based technique, for comprehensive characterization of the structural proper- 
ties of individual grains deeply embedded in polycrystalline samples. A de - 



taUed description of the technique is found in the book by H. F. Poulsen 1 2004ll , 

I2OOII 



and a detailed description of the geometry is found in iLauridsen et q;.I.I2001I1 . 



The method is based on the rotation method as described in section 12.4.11 but 
generally applied to polycrystalline samples. The sample is illuminated by 
a monochromatic beam, and the diffracted signal (in transmission mode) is 
recorded by use of an area detector behind the sample (as on figure IZ3l . Over- 
lap between reflections from different grains are avoided by matching the size 
of the angular rotations and beam size, with respect to the structural character- 
istics of the sample (grain size, sample thickness, mosaic spread and texture). 
Each detector image obtained includes diffraction spots from many grains, but 
each spot is generally only associated with one grain. 

The dedicated 3DXRD setup at the ID-11 beam-line at the European S5mchrotron 
Radiation Facility (ESRF, Grenoble) can be used in different modes: 

Spatial information on the shape of individual grains c an be prioritized for 
e.g. mapping of gr ain shapes during grain growth pPoulsen et ali I2OOII: 
ISchmidt et aqliool . 



Time resolution can be p rioritized for e.g. in-situ studies of grain growth 
llNielsenet a/.L 120031 1. 



Angular resolutio n can be prioritized fo r e.g. analysis of grain rotation during 



deformation iMargulies et a/.L 12001 1. 



Many grains in a polycrystalline sample can be investigated simultaneously 
due to the topographic approach to diffraction. The technique furthermore has 
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the advantage that only one (angular) degree of freedom has to be scanned 
during data acquisition, in contrast to spatial scanning techniques. 

The variant which is of relevance for the present study is what might be termed 
"far-field 3DXRD" where the detector is positioned so far from the sample that 
multiple full reflection families (Debye-Scherrer rings) are investigated. The 
orientation of the individual grains can be determined from such data. 

A fundamental part of the technique is the GRAINDEX program, by which it 
is possible to assign the individual diffraction spots obtained from a polycrys- 
taUine sample to individual grains. 

The individual reflections are first assigned to hkl families, and then by means 



of an algorithm known as "the grain digger" assigned to grains |Lauridsen et at 



200111. The output of the algorithm is a list of grains. For each grain the algo- 
rithm gives the orientation matrix, associated reflections, and statistical infor- 
mation on how confident the finding of the grain is. 

The most important statistical information are the uniqueness and completeness 
of the grain. The completeness measure how many of the theoretically acces- 
sible reflections have been found for the given grain. The uniqueness measure 
how many of the reflections assigned to a given grain that are also assigned to 
a different grain. 



2.5.2 The 3D crystal microscope 

The 3D crystal microscope is a white-beam (Laue) s5mchrotron-based tech- 
nique. This method is, as opposed to the 3DXRD method, a point-by-point 
measuring method using a micro meter si zed beam and what is known as differential- 
aperture X-ray microscopy [L arson et fl/.l,[2002 1. 

The dedicated setup at APS is described in details in lice et al. . 2005ll , and a 



recent review of the technique and results are given in Mice et al. 



2Q06I1. 



By use of X-ray optics the white-beam is focused to a sub micrometer spot 
(0.5 X O.Sixm). By scanning such a beam over the sample the structure can be 
probed in a 2D grid, but still probing all the material along the beam. The 
resolution in the depth is obtained by scanning an absorbing wire through the 
diffracted beam. Figure |23] illustrates the basic idea of the technique. By trian- 
gulation it is possible to reconstruct the Laue pattern from individual 3D voxels 
in the sample (the resolution in the direction along the beam is 0.5 - 0.7|^m). 

Local phase, crystal orientation and distortion of the unit cell can among other 
material parameters be gathered from the reconstructed Laue patterns. By 
modelling and simulation information can furthermore be obtained on the dis- 
location distribution. 

The major limitation to the technique is that the point-by-point data acquisition 
is very time consuming, and it can also be debated to what extent true bulk 
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Figure 2.5: Schematic of the 3D X-ray crystal microscope. The wire is moved in 
steps parallel to the sample to obtain spatial information along the beam. 



me asurements m ay be performed as the energy is rather low (results presented 
in I Ice et al. , 20061 go to a depth of 30|j.m in a Cu sample, but it is not explicitly 
stated that this is the depth limit). 

The technique has recently been extended to wh at is termed "scann ing monochro 



mafic differential-aperture X-ray microscopy" I Levine et al. , 2006|] . The tech- 
nique is specialized toward spatial measurements of elastic strains. It is, as 
indicated, a monochromatic technique. A tunable monochromatic beam (X- 
ray energy ~ 14keV) is used instead of a white beam, hence directly giving the 
lattice spacing in the illuminated sample part. The strain resolution is reported 
to be Aq/q = l-l-lQ-'^. 



2.5.3 3DXRD peak shape analysis 

The 3DXRD microscopy method has been generalized such that reciprocal space 
maps can be obtained for the individual reflections from deeply-embedded 
individual grains in a polycrystalline sample (a techniques which here will 
be termed "3pXRD peak shape analysis"). The technique was introduced in 
llPantleonet flZ.Ll2004ll . 

A specialized setup was developed at the 1-ID beam line at APS. By using of a 
narrow bandwidth monochromator and positioning the detector far from the 
sample a high resolution image was obtained in reciprocal space. By integrat- 
ing the stress rig into the setup it is possible to do in-situ measurements. 

The detector used was a 2D CCD detector, and by acquiring data while rotat- 
ing the sample over equidistant intervals (rocking) a 3D reciprocal space map 
could be generated of a reflection. As is general for the 3DXRD method the 
properties of the sample and beam were matched, so that the reflections from 
the individual grains did not overlap. 

Multiple reflections from the same grain can be found by use of the GRAIN- 
DEX program. The data acquired for this were taken using a large area CCD 
detector close to the sample (in the "far field 3DXRD" geometry). 
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In llPantleon et all 120041! measurements on 20 reflections from the same grain 
were characterized at tensile strains of 0%, 1%, 2.5% and 4.5%. The radial peak 
profiles (integrating in the azimuthal directions) show asymmetries, which in 
some cases are in line with what should be expected from the composite model 
(see page [191, but in other cases derivate from the expected. This is rational- 
ized in term of an anisotropy in the dislocation structure. The azimuthal pro- 
files were used for o btaining the orientation distribution function for the grain 



I Poulsen et fl/.L 1200511 . The method was further used for evaluatin g the effect of 



grain interactions during plastic deformation I Lienert et fl/.Ll2003l . 

Compared to traditional line broadening studies the technique has multiple 
advantages. The measurements are specific to the individual grains, multiple 
reflections can easily be investigated for the same grains and the data obtained 
are full 3D reciprocal space maps of the reflections. Compared to traditional re- 
ciprocal space mapping, the technique has the advantage that only one angular 
degree of freedom has to be scanned, due to the use of a 2D detector 

One of the ideas behind developing the setup was to gather in situ information 
on development of deformation structures. However, the interpretation of the 
reciprocal space maps still required heavy use of models. 

Just as the 3DXRD method avoids the averaging over multiple grains, it turned 
out that by enhancing the resolution in reciprocal space, this method could 
be turned into a direct probe for the individual subgrarns. This lead to the 
development of the "High Angular Resolution 3DXRD" method as described 
in the major part of this thesis. 



Chapter 3 



High Angular Resolution 
3DXRD 



The novel technique "High Angular Resolution 3DXRD" is presented in this 
chapter. 

The technique will be briefly described in section 13.11 including examples of 
the raw data and a presentation of the interpretation. The experimental setup 
is described in section 13721 followed by details on data taking and coordinate 
systems in section 13.31 and 13.41 The instrumental resolution is discussed theo- 
retically and experimentally in section |331 The different methods for analyzing 
the data are presented in section 13.61 The reproducibility of the results is dis- 
cussed in section ISTI followed by detailed argumentation for the interpretation 
of the data in section l3l8l High Angular Resolution 3DXRD is finally compared 
to other relevant techniques in section |3^ 

3.1 Overview 

This section gives an overview of the technique developed, the raw data ob- 
tained, the interpretation of the measurements, and a discussion of why a new 
technique is needed. The choice of samples is also discussed. This is a brief 
introduction and all issues will be discussed in detail in the following sections. 

3.1.1 The technique 

The aim of the technique developed is, as described in the introduction, to 
obtain high resolution reciprocal space maps of the broadened reflections from 
a deformed metal. 
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The experiments were performed at the 1-ID beam line at the Advanced Pho- 
ton Source (APS), where a unique combination of X-ray optics, detectors and 
mechanical setup enables 3D reciprocal space mapping with a high resolution, 
and a reasonable acquisition time. 

The meth od is developed on t he basis of the "3DXRD peak shape analysis" 
technique jPantleon et al. , 2003l described in the previous section (section |2. 5.31 1. 



The present setup allows for substantially higher resolution in full 3D in recip- 
rocal space, whereas the old setup was focused on obtaining low resolution 
peak shapes for a more traditional model-based analysis. 

The main difference between the setups, is that the present setup includes X- 
ray optics for focusing the beam onto a small area. This is needed if high res- 
olution maps are to be obtained in a feasible time. Furthermore the energy is 
different, and a different set of monochromators were introduced. 



The basic setup is shown in figure l3Jl The sample is illuminated by a monochro- 
matic X-ray beam, and the diffracted signal recorded on one of the two avail- 
able detectors. The sample is mounted in an Euler cradle allowing for rotation 
of the sample. This cradle is also used for rotation during exposure. 

Detector A is used for gathering "far-field 3DXRD" data (see section |2.5.1t . 
Such data are used either for finding the full orientation of diffracting grains 
by use of the GRAINDEX program, or for a manual search for interesting re- 
flections (see section l331 . 



Detector B is used for obtaining the high resolution 3D reciprocal space maps. 
This is done by rotating the sample around the x-axis in small consecutive in- 
tervals while obtaining data (the rocking method as described in section |2X3). 



3.1. OVERVIEW 
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Figure 3.2: Raw data. Each figure represents the intensity as recorded on the 
detector while rocking 0.007° around the x-axis. Each image covers a range 
of 0.6° horizontally (in t]) and 0.05° vertically (in 26). Read from left to right, 
top to bottom, the figures show consecutive intervals. The reciprocal 3D space 
map is generated by stacking such images. From iPaperlS 



Figure 13.21 shows 9 such consecutive images as obtained on the detector The 
sample is a Cu sample deformed to 2% in tension, and the investigated reflec- 
tion is a 400 reflection. A 3D reciprocal space map is constructed by stacking 
such 2D images. 

3D reciprocal sp ace mapping i s, as mentioned in section [2.4. 3i no new tech- 
nique (see e.g. I Fewster , 1997ll ). The advantage of the present setup is that 



due to the 2D detector a full 2D part of the reciprocal space is mapped per im- 
age obtained, in contrast to traditional pornt-by-point acquisition. The present 
method further has the advantage that the resolution is about equal in all three 
reciprocal space directions, and that the resolution can be changed by changing 
the rocking interval size and sample-to-detector distance. 



3.1.2 Interpretation of data 

From figure l3^ it is observed that the broadened reflection comprise a cloud of 
enhanced intensity upon which bright peaks are superimposed. The observ- 
able peaks are clearly separated in all three reciprocal space directions. 

Our interpretation of these structures, in the reciprocal space intensity distri- 
bution, is that the individual peaks arise from individual dislocation free sub- 
grains in the dislocation structure, and that the cloud stems from the dislocation- 
filled walls (the arguments for this interpretation are given in section |3l8l l. 

The radial position of the peaks are directly related to the mean strain in the 
scattering region, the width of the peak to the strain distribution within the 
scattering region, and the integrated intensity to the volume of the scattering 
region. 
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The techniques hence allows for a direct, model free, non-destructive investiga- 
tion of these properties, (mean strain, internal strain distribution and volume), 
of the individual subgrains. The properties can be investigated in-situ from 
dislocation structures deeply imbedded in individual bulk grains in a poly- 
crystalline sample. 

3.1.3 Samples 

The technique described in this chapter can naturally be used for obtarnrng 
data on single crystals but multiple issues (both experimentally and scientific) 
might suggest the investigation of polycrystalline samples. 

The use of polycrystalline samples has the advantage that the grain size will 
limit the investigated volume in the direction of the beam. If this volume is 
to large, the possibility of overlap between the peaks from the individual sub- 
grains becomes large, and it is impossible to separate them. In the case of single 
crystal samples the thickness would therefore have to be small, or the volume 
wo uld have to be restri cted by other means (a possibi lity is the use of a c onical 
slit llNielsenef g/.Ll2000ll or wire scanning techniques llLarson et aZ.1, 1200211 ). The 
use of polycrystalline samples has the further advantage that multiple grains 
with different orientations can be investigated under the exact same macro- 
scopic stress /strain conditions. 

The chosen material for most of the experiments was, as mentioned in the in- 
troduction, polycrystalline copper. The material used is 99.99% pure OFHC 
copper. The precise details of the samples vary between the individual exper- 
iments (see table WA) . Generally the material was cold rolled to a reduction 
of 80% to a final thickens of 300|j.m and then fully recrystallized by annealing. 
The material was characterized by electron backscattering diffractions (EBSD), 
and the mean grain size was foimd to be ~ 30 - 36ixm, resulting in about ten 
grains across the thickness of the sample. 

3.2 The setup at 1-ID (APS) 

The setup is illustrated on four figures: figure |3?T] which gives an overview of 
detectors, sample and stress rig, figure 13.31 which illustrates the X-ray optics, 
figure 13.51 (page |36l l showing a photograph of the main part of the setup, and 
figure l3l7l (page l4T1l which focuses on scattering geometry for reciprocal space 
mapping. 

The following main components will be described in detail below: 

• X-ray optics providing a focused beam with low beam divergence and 
low relative energy spread. Included is also a slit system defining the 
final beam size and beam position with respect to the sample. 
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Figure 3.3: Sketch of the X-ray optics setup. From Paper j. Supporting Online 
Material. 



A custom made stress rig, allowing for in-situ deformation. 

Translations and a Huber Euler cradle, allowing for positioning and ro- 
tating the sample. 

An area detector close to the sample for gaining "far-field 3DXRD" data 
for an overview of available reflections, also used for obtaining data for 
analyzing with the GRAINDEX program. This detector is mounted on 
a horizontal translation stage which allows it to be translated out of the 
diffracted beam. 

An area detector placed far from the sample, mounted on a vertical trans- 
lation stage, used for obtaining the high resolution reciprocal space maps. 



3.2.1 Optics and beam monitoring 

An overview of the optics developed by the A PS sector - 1 col laborators is p ro- 



An overview or tne optics developed by tne A rb sector - 1 col laborators is p ro- 
vided in figure |33I The optics consists of (see fPaper t\ and fS hastril, 120041): a 



pre-monochromator (two bent Laue crystals) [Shastri et ah, 20031/ & collimating 
refractive lens, a narrow-bandwidth mono chromator (two channe l cut crystals) 
and a set of saw tooth lenses for focusing ICederstrom et al. , 2002 1. 



By this combination of optics the beam obtained a unique combination of prop- 
erties: High energy (52keV), low vertical d ivergenc e (17|a.rad), small relative 



energy spread ^ = 7 • 10 ^, and high flux llPaper Jl . The horizontal beam di- 



T 

vergence is given by the source size and distance from this to the sample as no 
focusing exists in this directions [ Lienert , 2006], and is ~ 13|^rad. 

The low divergence and low energy spread allows for a high resolution of the 
reciprocal space mappings, and the high flux (due to the focussing) for a rea- 
sonable acquisition time. 

The final beam size impinging on the sample is defined by a set of slits posi- 
tioned close to the sample (~ 30cm). The intensity of the beam impinging on 
the sample is monitored by an ion chamber positioned behind (downstream 
of) the slits (seen on figure [ 
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Figure 3.4: Beam profile at the sample position in the horizontal (x) and vertical 
(y) directions. The profiles are acquired by scanning a micrometer sized LaBg 
grain through the beam and recording the scattered intensity. 

The position of the focal point of the focusing optics can be varied. The optimal 
flux would be obtained by focusing directly on the sample position, but as this 
in practice is difficult (since no beam monitor exists at this position) the focus 
point is either positioned at the slits before the sample, or on a high resolution 
detector which is positioned ~ Im behind the sample. In the later case, the 
result is a vertical focal size of the beam at the sample position of w 25|j.m. 
The total aperture of the focusing lenses is about 400|a.m and the integrated 
transmission w 80% [Lienert, 2006]. Taking into account the non-box-shape of 
the beam, the gain in total intensity is approximately a factor of 10. Using the 
slits the beam can be reduced to the vertical size needed. 



Measurement of the beam profile 

The beam at the sample position can be characterized by scanning a diffracting 
object through the beam and recording the scattered intensity. An example of 
a beam profile determined in this way is shown in figure 13.41 The profile is 
found by scanning a LaBg grain through the beam, the pre cise size of th e grain 
is unknown, but the mean size of the grains is w 2ixm I NISfl- bOOOl hence 
it is much smaller than the width of the beam. The slit size was nominal w 
20|J.m X 20ixm (approximately as the experimental run from which the data 
shown in figure l334l originate, suffered from problems with the motor drive of 
the slit blades). It is seen that the beam is box shaped in the a;-direction, and 
bell shaped in the y-direction; as expected. 

Beam stability 

The slits are fixed with respect to the diffractometer. This has the consequence 
that if the beam moves in the vertical direction the illuminated part of the sam- 
ple remains the same, but the total intensity and the intensity profile of the 
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beam at the sample position might change. Due to the long (w 21m) distance 
from the optics to the sample, even very small changes in angles will give de- 
tectable changes in vertical beam position at the slits, resulting in a change of 
the beam profile at the sample. Such instabilities have been a problem through- 
out this study. 

The changes in absolute intensity can easily be corrected for, by recording the 
integrated incoming intensity during data acquisition by the ion chamber be- 
fore the sampl^ More challenging is the change in beam profile, as this will 
lead to changes in relative intensities arising from different parts of the sample. 
A solution might be a feedback system on the optics, reducing the drift of the 
beam position. 

3.2.2 Stress rig 

The stress rig (shown in figure 13.51 is custom built to fit in the Euler cradle, so 
that in-situ experiments can be performed. The rig is displacement-controlled, 
and the deformation speed can be chosen in the range from very slow (ex- 
periments have been done as slow as ~ 4 • 10~^mm/s ) to rather fast (fastest 
experiment was done at w 0.6mm/s). 

The load on the sample is monitored by a load cell mounted below the sample. 
The load cell is read out by a control box, and readings are recorded by a log- 
ging system. The tensile stress on the sample is proportional to the measured 
load, as the gauge area of the sample is assumed to be constant. 

The strain is monitored by a strain gauge glued directly onto the gauge area 
of the sample. The strain gauge is read out using a Wheatstone bridge, and 
recorded by a logging system. 

3.2.3 Euler cradle 

The Huber Euler cradle allows for rotations around three axes. A 3-axis trans- 
lation is mounted on the inner rotation allowing for translation of the stress 
rig/ sample. 

The order of the rotations is such that (see figure l3^ on page|4T]for a definition 
of the rotations) cp is the inner rotation, x the middle rotation and lo is the outer 
rotation. 

On figure |3]5] the cradle is shown with the stress rig mounted. A cable guiding 
system is seen at the top of the image; this allows for free rotation of the stress 
rig, without human intervention. 

^The intensity logging system has been developed during this project to give the exact inte- 
grated intensity. In some of the early experiments different technical problems existed with this 
logging procedure, hence reported results do not in all cases have a perfect correction for intensity 
fluctuations. 
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Figure 3.5: Photograph of the main optical table showing: The stress rig 
mounted on the Huber Euler cradle, the Bruker CCD detector and visible 
through the Euler cradle the ion chamber. The insert shows a sample (this 
is an Al sample of different size than the Cu samples normally used) mounted 
in the stress rig. A strain gauge is visible on the sample. The load cell is the 
round device below the sample. 

Before measurements, the center of rotation of the cradle is made to coincide 
with the center of the beam. The centering is based on scanning a pin through 
the beam, and the geometry of the setup. 

However, for mechanical reasons the Euler cradle, with the mounted stress rig, 
does not have a totally well-defined center of rotation (a problem known as 
the sphere of confusion). This has the unfortunate consequence, that if a grain 
is centered in the beam at one (p, co, x setting, it might be off center at another 
angular position. The sphere of confusion can be measured in different ways, 
and corrections applied, however for most of the experiments presented the 
problem is minor, as the sample normally is only rotated a few degrees in the 
critical angles. 

The cv rotation (which always is aroimd the ac-direction), is used for the rocking 
procedure for gathering the high resolution 3D maps. An investigation showed 
that the drive could produce steps as small as 0.0005°, or better. 

3.2.4 Detectors 

The setup comprises two detectors (as shown on figure ISTTt . Detector A is used 
for obtaining low angular resolution data comprising many reflections from 
many grains (far field 3DXRD data). Detector B is the main detector used for 
obtaining reciprocal space maps. Detector A is moved out of the diffracted 
beam when data are to be acquired using detector B. 
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Detector A 

Detector A is a Bruker SMART 6500 CCD area detector. It is positioned ~ 0.3m 
from the sample. The exact distance varied from experiment to experiment, 
and was calibrated together with the position of the direct beam by a measure- 
ment on a known powder sample. The detector is close enough to the sample to 
allow the full 400 reflection family to be recorded when using copper samples 
and an X-ray energy of 52keV. The det ector ha s a pixel size of 161 \xm x 161 (a.m 



when used in 1024 x 1024 pixel mode llBrukerll . 



Detector B 

Detector B is a MarCCD 165 area detector, positioned at a horizontal distance of 
fa 3.9m from the sample. The exact position varied between the experiments 
and was found by calibration measurements. The detector is mounted on a 
vertical translation stage allowing for mapping different hkl's. 

The detector has been used in either 2048 x 2048 pixel or 1024 x 1024 pixel 
mode, the latter used for reduced readout time. The pixel size is 80.5i.im x 
80.5 (a.m and 161|.im x 161 (a.m respectively. The full width at half maximum of 
the point spread function is 100 |im LMARII . 



3.3 Selecting grains and reflections 

Before data can be obtained a reflection from an interesting grain (in the case 
of polycrystalline samples) has to be selected and centered on the detector. 

The grain is typical centered within the X-ray beam by spatial scanning of the 
grain with respect to the beam, and recording the diffracted signal. 

Two strategies for selecting reflections / grains have been employed in this study, 
and will be described below. The first strategy is based on a choice of reflection 
without full information on the orientation of the scattering grain. The second 
strategy is based on full information about the grain orientation as obtained 
from the 3DXRD technique (section [2.5.H . 



3.3.1 Reflection-based selection of a grain 

An overview over available reflections can be obtained by using detector A 
(close to the sample), and then simply scanning the orientation and position of 
the sample. By means of this a suitable reflection can be identified, typically 
based on a criterion on the angular separation between the chosen reflection 
and neighboring reflections. 
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Figure 3.6: Inverse pole figure shoving possible orientations of the tensitle axis 
in the crystallographic-coordinate-system, for known angles between the ten- 
sile axis and the {100} lattice plane normal. The angles shown are 0° (*), 2° 
(blue line), 7° (green line), 25° (cyan line), 45° (magenta line), 75° (yellow line), 
and 90° (red line). 

Based on one such reflection the corresponding grain can easily be centered in 
the beam in the x and y directions. In principle it is furthermore possible to 
clarify if the grain is within bulk or close to the surface. Several trigonometry- 
based methods have been devised but due to experimental uncertainties (mainly 
the sphere of confusion problem discussed in section |3.2.3t no perfect and sta- 
ble method exists (see section |3.3.3t . 

Based on one reflection it is possible to calculate the angle between the tensile 
axis and the normal vector corresponding to the scattering lattice planes. This 
leaves one rotational degree of freedom for the full orientation of the tensile 
axis with respect to the crystallographic coordinate system of the grain. In fig- 
ure l3.6| the possible orientations of the tensile axis, as represented in an inverse 
pole figure, are shown for multiple cases of 100 reflections with a fixed angle 
between the tensile axis and lattice normal. It is seen that in the case where 
the angle between the tensile axis and the {100} lattice plane normal is 0° full 
information about the orientation of the tensile axis is obtained, whereas in 
the case of 90° the orientation can be anjrwhere between (100) and (110). The 
orientation of the tensile axis is known to control the morphology of the mi- 
crostructure [Huang and Hansen, 1997] (see section[LT2J. 

This method for finding relevant reflections has been widely used in the present 
study as most investigations have been made on grains with the tensile axis 
close to a (100) crystallographic direction. 



3.3. SELECTING GRAINS AND REFLECTIONS 



39 



3.3.2 Grain-based selection of reflections 

Alternatively to working directly on reflections, a list of grains from the GRAIN- 
DEX program (as described in section l2.5.1t can be utilized. 

By working with the list of grains, and the list of reflections found, it is possi- 
ble to find grains and corresponding reflections which have certain properties. 
Generally the search is for grains of a certain orientation having interesting 
reflections, which are angularly separated from reflections of other grains. A 
number of MATLAB routines have been developed for this. 

When the full orientation is known it is furthermore possible to investigate 
multiple reflections from the sa me grain. This is very important as some exist- 
ing theoretical predi ctions (e.g. | Mughrabi/l 983*1) and experimental results on 
single crystals (e.g. lUngar et al. . 1984a]) show that differences exists between 



different reflections from the same crystal (see section l2.3.1t . 



3.3.3 Ensuring that a grain is in the bulk 

GRAINDEX-based methods 

The information from the GRAINDEX program can be used in two ways to 
find individual grains which are in the bulk of the sample. 

The first method is based on the observation that the completeness parameter 
associated with each grain found (see section [2. 5. Il l is unlikely to be high for a 
grain not in the bulk of the sample. This is true under the following conditions: 

• The sample is positioned such that the center of the sample is at the center 
of rotation of the Euler cradle. 

• The vertical beam size used is much smaller than the thickness of the 
sample. 

• The grains are substantially smaller then the thickens of the sample. 

• The dataset obtained for the GRAINDEX analysis covers a substantial 
angular range. 

A grain not in the bulk of the sample will rotate out of the beam at some point 
during the data acquisition, if these requirements are met. The consequence 
is that a number of expected reflection will be missing, and GRAINDEX will 
therefore return a low completeness. 

Secondly, the knowledge of the full orientation of the grain, allows for selecting 
multiple reflections from the same grain, which are separated by some rotation 
around the y-axis. The full position in the {x, z) plane can be found by scanning 
the grain through the beam at two such different orientations. It can therefore 
be directly determined if the grain is in the bulk of the sample. 
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Single reflection-based methods 

If only one reflection from a grain is known it is still possible to find different 
angular settings of the Euler cradle which gives access to the reflection. It is 
therefore possible by triangulation to find the full spatial position of the grain, 
and hence to determine if it is in the bulk. 

A simple example is that a grain rotated by an angle of 29 around an axis per- 
pendicular to the scattering plane will bring "minus" the scattering vector into 
scattering condition. 

The problem of such methods is that the sphere of confusion makes it difficult 
to determine if the grain is in the bulk of the sample. 

3.4 Reciprocal space mapping 

The general considerations about obtaining 3D reciprocal space maps will be 
presented in this section and a useful reciprocal space coordinate system will 
be defined. 

Figure l3l7l is general for this section. It defines the real space coordinate sys- 
tem, available rotations angles and illustrates the reciprocal space coordinate 
system. 

3.4.1 Data collection 

As shown on figure I3l7l the sample is placed in a 3-axis rotation stage and illu- 
minated by a monochromatic X-ray beam. The diffracted signal (in transmis- 
sion mode) is recorded on an area detector positioned on a vertical translation 
stage behind the sample. 

Before data can be obtained from a desired reflection, from a specific grain, 
the sample is rotated, and the detector translated, such that the reflection gives 
rise to a diffraction signal on the detector. That is, the sample is rotated such 
that the diffraction condition, for the relevant reflection, is fulfilled with the 
scattering vector close to the {y, z)-plane, and the detector translated to the 
relevant height for the scattering angle of the reflection. 

Obtaining reciprocal space maps 

Data are (as described in section l3.1.1t obtained by rocking. By this a 3D part of 
reciprocal space is mapped onto the detector, by integrating out one dimension. 
In the present setup the co rotation is used for the rocking (as it is approximately 
perpendicular to the scattering planes for the investigated q- vectors). The 
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Figure 3.7: Sketch of the experimental setup. Defining the laboratory coordi- 
nate system x, y and z, available rotations of the Eulerian cradle (p, x and to 
(with ^ the inner rotation and co the outer), angles 20 and ?/, and reciprocal 
space coordinate system qx, qy and qz- In the indication of the angle between 
the (/z-axis and the horizontal plane, 9q is the scattering angle of the unde- 
formed sample, and co is to be read as the present value of this rotation with 
respect to vertical. L is the horizontal distance between the detector and the 
point of rotation of the Eulerian cradle. The tensile axis is along the y-direction 
for CO = Q, and x = ^- From LPaper ILL 
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rocking interval is termed Aco and one image on the detector is termed an o;- 
slice. By acquiring and stacking adjacent a;-slices a 3D map of a part of the 
reciprocal space is obtained. Each point in this dataset represents the integrated 
intensity of a voxel in reciprocal space. 

The rocking interval can be chosen freely (within the physical limits of the 
setup) and hence the resolution of the 3D map in one azimuthal direction 
changed. The resolution in the other directions can, in principle, be changed 
by changing the sample-to-detector distance, but this option was not used in 
the present investigations. Data are generally acquired with three different res- 
olutions in the one azimuthal direction: 

• High resolution; where the rocking interval size is matched to the other 
azimuthal voxel size, as given by the pixel size on the detector and sample- 
to-detector distance (see section 1331) 

• Medium resolution; Resolution matched to give ~ 10 w-slices over the 
azimuthal spread of the reflection investigated. 

• Low resolution; integration over the full spread of a reflection, mainly 
used for obtaining radial intensity profiles integrated over the full az- 
imuthal spread. 

3.4.2 Definition of the reciprocal space coordinate system 

In this secti on the sca ttering geometry is described in detail (parts of this is also 
included in IPaperll) , giving the relevant connections between the real space 
coordinate system, angles, and reciprocal space. The notations are generally 
defined on figure |371 

The connection between the positions in the real space coordinate system, as 
given by x, y, z, and the angles 29 and rj are given by: 



q in the laboratory system 

The direction of the incoming beam is by definition along the positive z-axis, 
hence: 



= arctan (— x/y) , 



(3.1a) 



26 




(3.1b) 



fco = ^(0,0,1). 



(3.2) 
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The outgoing beam is given as: 

f-sm{rj)sm{2e) 
s(7/) sin(~'" 
cos (20) 



fc = — I cos(7/)sin(20) I . (3.3) 



Combining these two equations the scattering vector (q-vector) in the labora- 
tory system is given as: 

2^ /-sin(/7)sin(20)^ 
'/lab = -Y- cos(;7)sin(20) | . (3.4) 
\ cos(20)-l 



q in the q-space coordinate system 

The q-space coordinate system (see figure l3?7t is defined with respect to a per- 
fect crystal having the scattering conditions fulfilled at current x and (p settings 
wittQ a; = 0, and the scattering vector in the {y, z)-plane. This leaves one an- 
gular degree of freedom for the coordinate system. This is fixed by aligning the 
<jx-axis with the real space x-axis. This has the advantage that qx then becomes 
invariant under rotation in lo, which simplifies the transformation equations 
significantly. 

This definition further has the advantage that small angle approximations can 
be used for two reasons. Firstly the scattering angle for a deformed crystal, 
6, is close to the scattering angle for a perfect crystal, 6q, as the elastic strain 
is small. Secondly, the scattering vectors for the reflection investigated are all 
close to the (y, z) -plane, having small rj, because of the geometry of the setup. 

The relation between the <j -vector represented in the laboratory system, qi^^,, 
and the representation in the q-space coordinate system, q, is given by: 

^10 ° \ 

cos{9o + co) - srn(0o + <7iab- (3-5) 
yO sm{dQ + Lo) cos{6q + co) J 

Combined with equation 13.41 this gives the following expression for the q- 
vector in the reciprocal space coordinate system, as a function of scattering 
angles: 

— sm{i]) sin(20) 
cos(0o + ^) cos{r]) sin(20) — sin(0o + cos(20) + sin(0o + ^) 
^sm{6o + co) cos{rj) sin(20) + cos(0o + ^) cos(20) — cos(0o + <^)/ 

(3.6) 



^If cv is read as relative to some arbitrary zero point, this condition is relaxed. However, in 
the present studies, focus is on cases having the tensile axis close to the scattering vector, and it is 
therefore feasible to keep the definition of cv absolute. 
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Assuming ?/ w 0, a; w and w equation 13.61 can be expanded to first order 
in the small quantities to: 



where \qo\ = ^sin{eo). 

This equation is the one used for all calculations in this thesis, as the small 
angle approximation is appropriate in all cases. 

3.5 Instrumental resolution 

The true shape of a reflection arising from the material properties, are distorted 
by the intrinsic instrumental resolution function. The instrumental resolution 
is, in this section, discussed on the basis of the parameters of the setup, and on 
the basis of measurements. 

3.5.1 Relations between experimental parameters and resolu- 



Equations for the contributions to the instrumental resolution function in re- 
ciprocal space are, in this section, presented as a function of parameters of the 
setup. 

Relation between detector pixel size and q-space resolution 

The connection between a position on the detector and the corresponding an- 
gular position is given by equation 13. II To connect a given size of a real-space 
region (e.g. a pixel) to a size measurement in reciprocal space an approximate 
equation is derived. The assumption is that the length scale of the region is 
small, and that the region has one corner at a reference point at x = 0, y = yo, 
z = Zq. The region is furthermore assumed to be in the (x, i/)-plane. The ref- 
erence y-position can be expressed in terms of a reference scattering angle as: 
tan(20o) = yo/zo- 

Expanding equation I3.1al for small x values, A^, and equation I3.1bl for small 
changes in y, Ay, gives the following approximations for changes in angle, A9 
and At], with changes in real space coordinates: 




(3.7) 



tion 



At] 



Ay 1 ^ Ay 

zo 2(l+tan2(20o)) ~ 2zo 
A.T 1 A.^ 1 



(3.8b) 



(3.8a) 



Zo tan(20o) 2zotan(0o)' 



where the second expansion is to first order in 6q. 
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Combining this with equation 13.71 gives the following approximation, to first 
order in 6q, for the dimensions in reciprocal space corresponding to a real space 
region with size A^, A^: 

A<?,,«J^^. (3.9a) 
ztant^o 

In the <jz direction the voxel dimension is determined by the rocking interval, 
Aco, from equation 13. 71 

Aqz = I (JO I Aw. (3.9b) 

These equations (|3.9a| and |3.9bt can be used to calculate the voxel size in q- 
space from the detector characteristics such as pixel size and point spread func- 
tion, combined with the rocking interval. 



Relation between beam properties and q-space resolution 

The beam is characterized by the relative energy spread AE/E, the vertical 
beam divergence and the horizontal beam divergence. The first two quantities 
mainly limit the resolution in the <jy-direction. The horizontal divergence is 
mainly related to the resolution in the qx direction. 

The relation between energy and wavelengthH gives: 

AA AE 

Ao to 

which combined with a Taylor expansion of Bragg's law, and the expression 
for qy (from equation |3.4t gives: 

Aflu AE 

T^ = ^' (3-11) 

which directly gives the contribution to the resolution function in the radial 
direction from the relative energy spread. 

The vertical beam divergence a. (normally expressed as full width at half max- 
imum) directly adds to the scattering angle 26, the spread in cjy is therefore 
foimd from equation | |3.7t as: 



^The wavelength is connected to the energy of the X-ray beam, E by: A — 'J^, conveniently 

12.398 
£[keV] 



expressed as A [A] = ■'^■^'^ 



46 



CHAPTER 3. HIGH ANGULAR RESOLUTION 3DXRD 



The horizontal beam divergence (expressed as the angle /3) gives a broadening 
on the detector in the horizontal direction of Ax = L • /5 (L being the horizontal 
sample-to-detector distance), which according to equation | |3.9a| | corresponds 
to a spread inqx of: 

2tan(t^o) 



3.5.2 Calculated instrumental resolution 

In this section the contributions to the instrumental resolution from the differ- 
ent parts of the setup are collected and reported in reciprocal space units. It is 
further discussed how these contributions should be added. 

The conversion is based on the equations from the last section, the experimen- 
tal parameters from section |372l a sample-to-detector distance of L = 3856mm, 
and a copper 400 reflection at an X-ray energy of 52keV. The corresponding 
scattering angle is lOg = 15.16°, and length of the reference reciprocal lattice 
vector is qo = 6.9525. 

The results are reported in table 13.11 It is seen that the resolution-defining fac- 
tors all give rise to similar instrumental broadening in reciprocal space. The 
resolution in the i^z-direction is given by the rocking angle, and can be matched 
to the resolution in the other directions. 

Three contributions limit the resolution in the radial direction: the point spread 
function, the vertical beam divergence, and the energy spread of the beam. It 
is however not clear how the different contributions to the resolution should 
be added. If the profiles are Lorentzian, full widths at half maximum should 
be added linearly and if Gaussian the contributions should be added quadrati- 
cally. As the real profile probably is between these two extrema we can provide 
a lower and upper limit for the theoretical experimental broadening in the t/y- 
direction as follow: 

A^j/^max = (7 + 4 + 5) -lO^U-i = 16- 10"^-^ (3.14a) 
A(jy,i„in = ^ (72 + 42 + 52) • lO-^A-i = 9.5 • lO-^A-i. (3.14b) 

In the <j^-direction the point spread function and horizontal beam divergence 
give the theoretical broadening: 

A^^,max = (7 + 3)-10"^A"^ = lO-lO-^'A"^ (3.14c) 
A(?;,,mm = ^ (72 + 32) . lO-^A-i = 7.6 • 10-^-1. (3.14d) 
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Given resolution 


Reciprocal space 


Detector pixel size|^ 


80.5i^m X 80.5(a.m 




= 5 


10 




-1 


Detector point spreac^ 
Relative energy spread^ 


100nm X lOO^m 




= 7 


10^ 




-1 


7 -10-5 


Aqy 


= 5 


10 


-4a- 


-1 


Vertical Beam divergenc^ 


17|j.rad 


My 


= 4 


10- 


-4A- 


^1 


Horizontal Beam divergence^ 


13|J.rad 


Aqx 


= 3 


10" 




-1 


Rocking angle£ 


0.0005° 


Aqz = 


0.6 


10- 


-4A- 


-1 



"Section lalZil [MAR]. 

''Section lam IPaperl. 

■^Section lsn] rLienertl l2006h . 

''Minimal rocking angle tested on the Huber Euler cradle. 



Table 3.1: List of nominal measures of resolutions, and the corresponding re- 
ciprocal space widths. Calculations are done under the assumption of a 52keV 
X-ray beam, a sample-to-detector length of L = 3856mm, and a vertical posi- 
tion of the detector corresponding to a 400 reflection from Cu. 



3.5.3 Measured instrumental resolution 

The following investigation does unfortunately not lead to a full characteriza- 
tion of the experimental broadening; the analysis will only indicate how close 
the experimental broadening is to the theoretical limits. 

The experimental broadening was investigated by measurements on single 
scattering domains in LaBg powder. The measurement was performed using 
a beam of w 10|j.m x 8ii.m (similar investigations using a larger beam indi- 
cates that the effect of the beam size is insignificant). Thirteen well separated 
peaks belonging to the 421 hkl family (w hich has a d-spacing close to that 
of the 400 planes of copper jMSlt I2OOOII ) was identified from a large recip- 



rocal space map (100 images using a rocking angle of 0.001° (corresponding 
to Aqz = 1.2 • lO^^A^^). For each peak, profiles along approximately the (\x, 
(]y and qz directions were obtained and fitted individually to a pseudo-Voigt 
function plus a constant term (see section l3.6.2l and l3.6.3l for discussion of this 
procedure and the function used). 

Figure l3!8l shows the resulting individual widths and the mean values. Also 
plotted is the theoretical size broadening fro m the LaB^ g rains Aijsize = 0-26 • lO^'^A- 
based on the mean size of the grains, 2 |j.m fNISl^, 120001, and the Scherrer con- 
stant for a spherical grain, 0.829 [Langford and Wilson, 1978]. It is observed 
that the broadening in the (]x and (jy-directions are very similar and substan- 
tially larger than the size broadening. The broadening in the ^z-direction is on 
the other hand very close to the size broadening. 

The mean values and differences, linearly and quadratic, from the theoretical 
size broadening are given in table 13.21 These numbers for qx and <jy can be 
compared to the limits given in equation 13.141 
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2 4 6 8 10 12 14 
peak nr 

Figure 3.8: Peak width of thirteen 421 reflections from individual scattering 
LaBg domains. The sample-to-detector distance was 3736.7mm, and a vertical 
detector position corresponding to the 400 reflection family of Cu at a X-ray 
energy of 52keV was used (vertical center position of the detector was 1095mm 
above the direct beam). Widths in the three directions of the individual peaks 
are plotted as: Aqx (x)/ (o) and Aqz (D)- Averages over the widths in the 
three directions are plotted as: {Aqx) (dashed line), {Aqy) (dash-dotted line) 
and {Aqz) (dotted Une). The theoretical mean size broadening (equal in all 
three directions) is plotted as: (Assize) (full line). 
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{Aq) 


(Aq) - {Aqsize) 

[10-4A-1] 


{Aq)2 - {Aqsize)^ 

[10-4A-1] 




11.3 


8.7 


11.0 


qy 


11.6 


9.0 


11.3 




3.7 


1.1 


2.6 



Table 3.2: Mean width of the 13 reflections, as shown in figure 13.81 in the 
three reciprocal space directions, {Aq). Included is the linear and quadratic 
differences between the measured widths and the theoretical size broadening 
(A^jsi^e) =2.6-10-4A-i. 

The range of the measured experimental broadening in the (/y-direction over- 
lap with the theoretical calculated range, and indicates that the real experimen- 
tal broadening is close to the theoretical lower limit. 

The range of measured experimental broadening in the t/x-direction likewise 
overlap with the theoretical range, but in this case the measurement indicates 
that the real experimental broadening is close to the upper theoretical limit. 

The observed experimental broadening in the (/z-direction is very close to the 
resolution as defined by the rocking angle interval. It is therefore not possible 
to evaluate the intrinsic broadening in this direction, but only to observe that 
it is seen to be small. 

The analysis above only set limits on the experimental broadening. If a true 
deconvolution is to be performed, a better method for characterizing the ex- 
perimental broadening has to be devised. In the following scientific cases it 
will however be shown that much can be learned form the simple estimates of 
the experimental broadening. 

3.6 Data analysis 

One or multiple 3D intensity distribution maps are obtained for each exper- 
iment. These datasets can be seen as 4D movies. The major challenge with 
respect to data analysis is how to gather relevant physical information from 
such large dimensional datasets. 

As the broadening of the reflection is mainly in the azimuthal directions, the 
major technique for getting an overview of the intensity distribution is a pro- 
jection onto the azimuthal plane. 

Projections in one azimuthal direction over small parts of the perpendicular 
azimuthal direction have also been used, in cases where the radial (strain) in- 
formation was analyzed. 

Such 2D projections give a good overview of time/strain dependencies when 
arranged in tables. 
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Projections onto the radial direction carry information on the distribution of 
strain in the illuminated part of the grain (equivalent to traditional line broad- 
ening measurements). Such projections have been used for defining mean 
strain in the grains, and for comparing to the present line broadening-based 
results. 

The background for the projections are rebinning of detector-pixel-data to the 
reciprocal space. The rebinning method used is presented in section |3.6.1| Sec- 
tion |3.6.2l presents the projection methods. 

The main results of the technique are the properties of the observed single indi- 
vidual peaks. The analysis method used on these peaks is presented in section 
13331 

To gather statistical information on the observed cloud of enhanced intensity, 
and further information on the radial position of the individual peaks, a statis- 
tical analysis method has been devised, this is presented in section |3.6.4| 

The general technique of converting the integrated intensity to volume of the 
scattering entity is finally presented in section |3.6.5l 

3.6.1 First order rebinning 

The problem of converting data obtained as pixels on a 2D detector to pixels in 
reciprocal space, or equivalent to angular space defined by pixels in a rj, 26 co- 
ordinate system, is general for many applications of area detectors. The major 
problem is that a square pixel on the detector does not corresponds to a square 
pixel in reciprocal space, and that the shape and size of the pixel in angular 
space depends on the position of the pixel on the detector (figure l3l9l illustrates 
the differences in the case of the present setup). 

Multiple software packages (such as |'FIT2D']) exist that are able to perform 
the general transformation, known as rebinning. The geometry of the present 
setup does, however, seems to give problems with FIT2D as the beam center 
is located far outside the image obtained. A first order rebinning scheme has 
therefore been used. 

The differences in area in reciprocal space between the different pixels are mi- 
nor, as the detector (and more important the interesting reflection) only covers 
a small angular range in the present case. This can also be seen on figure I3l9l 
where the differences are illustrated in the extreme case. 

For the typical case (as represented in the data presented in lPaperlB) the max- 
imal difference in the size of the voxels that corresponds to the pixels in the 
relevant region was calculated. In the azimuthal direction, qx, the change in 
length is at most 0.33% of the minimum length, and in the radial direction, qy, 
the change is at most 0.05% of the minimum length. In the other azimuthal 
direction q^ the size is constant. This gives a maximal difference in the volume 
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Figure 3.9: The mapping of one pixel squared on the detector onto the {qx, ciy)- 
plane for three horizontal pixel positions. Geometrical parameters are chosen 
within typical values of the actual setup. The mapped pixel are located at (left 
to right) pixel nr. —1000, 0, 1000 with respect to the horizontal beam center, 
and vertically close to the position of the 400 reflection for Cu. 



of the voxels of 0.4%, and a first order rebinning scheme, which does not take 
this volume difference into account, is therefore appropriate. 

The rebinning is not performed separately but is integrated into the projections 
as described in the following. 



3.6.2 Projections of reciprocal space 

Based on a data-set that have been rebinned to reciprocal space, projections are 
trivially done as a sum over voxels in the appropriate directions. 

The general "first order rebinning" based projections are performed in the fol- 
lowing way: 

• Define the bins of the final projected intensity distribution: 

The bins are defined using the average pixel position in the selected re- 
gion, in the x or y-directions for the projections onto the radial direction 
or azimuthal plane respectively. 

• For all voxels in the data set, calculate the q'-position. 

• For all voxels find the closest bin in the final projection, add the intensity 
to this bin, and record that a voxel has been added to the bin. 

• For each bin in the projection normalize by the niimber of added voxels 
to that bin. 



Azimuthal projection 

Projections onto the azimuthal plane have been done in two ways; the proper 
projection based on first order rebinning as described above, and an even sim- 
pler projection, based directly on the detector pixels. 
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The latter case is based on the fact that the horizontal plane of the detector is 
approximately in the (jt-direction, hence a good approximation to the projec- 
tion onto the azimuthal direction is obtained by summing over the columns of 
the detector pixels. 

In the (/.v-direction this is very accurate, as the change in (jv-position of the 
pixels over one column is much smaller than the bin size of the final projection. 
However, the <jz-position of a pixel depends directly on the 26 position of the 
pixel (according to equation |3.7| | and it is found that by just summing over 
columns some voxels might end up in the wrong bin in the projection. The 
error is however minor, and no significant differences are observed between 
the two projection t5^es. 

Generally the proper projection has been used for the large azimuthal projec- 
tions (as the one in figure 3 in Paper ID , and the simple projection for small 
projections (as the ones presented e.g. for reproducibility in section |3l7t . 

Radial projection 

All integrated radial peak profiles reported in this study (as the one in figure 
14. H have been calculated using the first order rebinning technique. 

The radial peak profiles are normally fitted (by least square fitting) using ei- 
ther a split pseudo-Voigt or normal pseudo-Voigt function. The pseudo-Voigt 
function is a weighted sum of a Gaussi an and Lorentzian function (for a discus- 
sion of this function see |Enzo and SchiffiniLri9 99]). The normal pseudo-Voigt 
function has one parameter describing the width of the peak (in the version 
used the parameter is the full width at half maximum). In contrast the split 
pseudo-Voigt function has two width parameters, one giving the half width 
of the lower (with respect to the maximum) part, Aj, and one giving the half 
width of the upper part, A2 (see figure l3?T0t . This can therefore be used for 
fitting asymmetric peaks. The two width parameters for the split pseudo-Voigt 
are normally added (A^ + A2) to give the full width at half maximum, and 
subtracted (A^ — A2) to give a measure of the asymmetry. 

As quantitative results are derived from the radial peak profiles, it is important 
to know if the first order rebinning scheme distorts the peak profile. 

To investigate this a detector image of a known radial peak profile was calcu- 
lated. The peak profile is assumed to be constant in the azimuthal directions 
and to be described by a split pseudo-Voigt function in the radial direction. 
Typical parameters for the geometry of the setup and peak profile were used. 
To simplify the calculation the parameters of the peak profile, were converted 
to angular units, and the peak described as function of the scattering angle, 
f{20). The intensity in each pixel, /pixel/ was calculated as: 




(3.15) 
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Figure 3.10: Illustration of the split pseudo-Voigt function. Indicated are the 
two "half width at half maximum" parameters ( Aj and A2), maximum position 
(xmax), maximum value (/max)/ and half maximum (/max/2). 



where the functions 29{x, y) and rj (x, y) are given by equation 13. II The integra- 
tion was carried out numerically. 

The simulated image was treated in the same way as the normal data, and the 
simulated parameters and parameters found by the fitting were compared: 



Parameter 


Relative difference between fit- 




ted and simulated parameter 


Peak position 


1.5 • 10~*% 


Full width at half maximum 


0.15% 


Asymmetry 


0.7% 



From this it is concluded that the first order rebinning scheme does not intro- 
duce any significant errors in the relevant parameters of the resulting radial 
peak profile. 



3.6.3 Single peak analysis 

Individual peaks have been identified either from the raw 3D data sets, or more 
frequently from azimuthal projections. 

The position of the maximum, q, the full width at half maximum (FWHM) in 
the three directions. At/, and the integrated intensity, /, have been found for the 
individual peaks. 
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Figure 3.11: Example of fitting a single peak to three individual pseudo-Voigt 
functions along the three reciprocal space directions. Dots indicate measured 
intensities and lines the f it obtain ed, the horizontal line indicates the full width 
at half maximum. From lPaper III 



The analysis is based on the full 3D intensity distribution in reciprocal space, 
according to the following procedure: 

• The voxel with maximum intensity was identified. 

• Three one-dimensional intensity profiles were gathered (along approxi- 
mateljQthe qx^ qy and (jz -directions) going through the maximum voxel, 
and not integrating in the other directions. 

• A one-dimensional pseudo Voigt function plus a constant term were fit- 
ted (using least square fitting) independently to each of these three pro- 
files (figure |3 . 1 1 1 show an example of such a set of three fits). 

The fits directly give the position of the maximum in reciprocal space and full 
width at half maximum of the peak in the three directions. 

On figure |3J11 it is seen that the constant level in the two azimuthal directions 
{qx and qz) are approximately the same, this indicates that the procedure in 
these directions correctly separates the cloud and the peak. It is also seen that 
the constant level in the radial direction {qy) is lower, in fact it is almost at the 
detector background level. The reason for this is that it is very hard to separate 
the cloud and the peak in this direction as both tend toward 0. This has the 
consequence that the method has a tendency of overestimating the width in 
the <7y-directions. The uncertainty on the width is mainly attributed to the 
uncertainty on the level of the cloud, and hence on the half maximum level. 
It can be estimated individually from the fits, by varying the constant level 
between the fitted values. 

The integrated intensity in a peak is found from the fits to the peaks. It is as- 
sumed that the average constant level and peak height found in the azimuthal 



'*The lines are taken along straight pixel lines on the detector. 
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directions can be used in the radial direction. The uncertainty of the total inten- 
sity was estimated by varying the background level used and the peak height 
between the minimum and maximum fitted values and constructing the corre- 
sponding upper and lower limits of the intensity. The estimated uncertainty is 
w 15%. 

3.6.4 Statistical analysis of intensity distribution 

By the methods presented in the previous sections | |3.6.2| and |3.6.3t information 
can be gathered on the overall strain distribution in the illuminated grain, and 
on the strain in the individual scattering entities giving rise to the peaks. In- 
formation on the strain distribution in the material which gives the cloud is a 
very good complement to these data on the grain and the subgrains. 

The problem with gathering information on the cloud is that it is very hard 
to locate parts of the cloud which have no influence from the superimposed 
peaks. A statistical method for analyzing the overall distribution of strains in 
the illuminated material has been devised to overcome this problem. 

The algorithm of the technique is as follows (see also figure |3?T2)| : 

• All possible radial^ line profiles are gathered (that is, one for each pixel /bin 
in the azimuthal plane). 

• For each profile the maximum intensity, Imax, is found. 

• The (/y-position of this maximum, cjy^max, is found. 

• The((;y, max/ ^max) pairs are bmned into log-intensity bins of constant width. 

• The mean tjy^max value is finally calculated for each bin, together with the 
number of points in each bin. 

An overview of the strain distribution is gathered by plotting the mean (jy max 
as a function of intensity (an example of such a plot is given in figure lXTSl . The 
points at low intensity will correspond to an average over mainly the cloud, as 
the cloud has a much lower intensity than the peaks. In this way information 
on the average strain in the material giving the cloud can be obtained. 

3.6.5 Volume calibration 

The integrated intensity from a scattering entity is, as described in section l2.4.1[ 
proportional to the volume of the scattering entity, under the assumption of 
kinematical scattering. This gives a possibility for finding the volume of both 



^To simplify the procedure the line profiles are gathered along detector pixel coliumis, which 
are all very close to the radial direction. 
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Figure 3.12: Ilustration of the statistical analysis method. The left part shows 
a raw detector image. The right part is the intensity profile along the red line 
indicated on the raw detector image. I,„ax and qy,max are indicated on the plot. 



the part of the grain illuminated and the volumes giving rise to the individual 
peaks. 

The proportionality constant is found by acquiring a diffraction pattern form 
a reference powder sa mple of the same material, using the same hkl family 
I Lauridsen erall,l2000|l . 



The general equation for the total diffracted intensity per angular unit from 
a powder was given in section 12.4.21 Equation 12.281 (on page |22) was derived 
under the assumption of an unpolarized X-ray beam, and a re-derivation is 
needed to generalize to other sources. 

In the present case two facts simplify this need for a modification; the scattering 
angle is small due to the high energy and the reference measurement is only 
carried out over a small part of the Debye-Scherrer ring close to f] = 0. Both 
facts have the consequence that the polarization factor is very close to the one of 
vertical scattering from a horizontal polarized beam, which is constant, P = 1 . 
Equation 12.281 can therefore be used with P = 1 as a model for the reference 
intensity measured on the detector. 

The normalized reference intensity, Epg^der' defined as: 

, _ £powder4sin(gM/) -360° 

'^powder ~ t t/ niM ' (oAb) 

^O,powder ^powder'"'-^' 

with /o,powder being the intensity during the reference measurement and other 
symbols as in equation |2.28l 

The volume of a scattering entity ( K^ystai) can then be found as: 

„ _ Ecrystaia^sin(20;^ 

''crystal ^ J~F/ W-^'') 

powder 

where Ecrystal is the integrated measured intensity of a reflection associated 
with the entity, Jq the incoming intensity during the measurement and other 
symbols as in equation l2.27l (on pagel22l. The term 1 / 1 sm{rj) \ in equation l2.27l 
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it replaced by 1, as the rotation axis is perpendicular to scattering plane. It is 
assumed that the same material has been used for the reference powder as the 
one investigated, and that the powder measurement is carried out for the same 
reflection family as the one investigated. 

In practice a highly deformed, cross rolled sample of the same metal is used 
as the reference material. A large slit opening is used for the measurement 
such that the diffraction pattern represents the average over sufficiently many 
grains that X-ray powder theory can be applied. Furthermore, the sample is 
put on a spinner to eliminate textural effects. 

The biggest problem with this volume calculation scheme is that the incoming 
beam intensity has to be known, and the beam is, as was discussed in section 
13.2.11 vertically bell shaped, giving a rather non-uniform incoming flux to dif- 
ferent parts of the illuminated volume. 



3.7 Reproducibility 

The reproducibility of the 3D intensity maps has been tested many times, by 
repeating the mapping of a specific area of reciprocal space. The stress rig does 
only operates in displacement control, some stress and strain relaxation might 
therefore be observed during the measurements. 

Figure l3.13l shows a reproducibility study conducted on a deformed Cu sam- 
ple in connection with Paper I (full details on the experimental conditions are 
given in the paper). A 3D map of the same region of reciprocal space was ob- 
tained three times at a tensile strain of 3.3%. Each map took ~ 9 minutes to 
obtain and was taken immediately after each other. 

The logging procedure for the incoming intensity was unfortimately not work- 
ing properly when this dataset was obtained. Therefore no intensity normal- 
ization has been applied; the total scattered intensity in the three maps are 
however almost constant (within 0.3%), indicating a stable beam. 

It is seen that the three maps comprise the same features, with only minor 
variations of the exact intensities. This demonstrates that the method gives re- 
producible results, within the time limits of these investigations, w 30 minutes. 

A longer investigation consisting of 40 repeated mappings, was conducted in 
connection with the study presented in Paper VI. The total time spanned was 
~ 1.5h, and data sets were obtained immediately after each other The datasets 
were normalized with the total incoming beam intensity, and the scattered in- 
tensity was constant within ±6%. During the experiment a small stress relax- 
ation was observed (the stress drop is w 1.4%, and strain increment w 0.2%). 

A set of azimuthal projections are shown in figure l3T4l representing every sixth 
map. It is seen that the general features of the intensity distribution are stable 
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Figure 3.13: Reproducibility test in connection with iPaper 1 . A 3D reciprocal 
space map was obtained three times of the same part of the reciprocal space 
near a 400 reflection from a Cu sample deformed to 3.3%. Each part of the 
figure is a projection onto the azimuthal plane. Each map took w 9 minutes to 
obtain and was taken immediately after each other. 



on the time scale of the experiment. Some minor changes are observed, spe- 
cially in the low intensity peaks, these effects are likely to arise from changes 
in the beam profile illuminatrng the sample, but could also be due to physical 
changes in the sample due to stress relaxation. 

These investigations demonstrate that the observed structures in the intensity 
distribution are reproducible. 



3.8 Detailed arguments for the interpretation 



The interpretation of the bright isolated spots as arising from individual dislo- 
cation depleted subgrains and of the cloud as arising from the dislocation-filled 
walls was introduced in section [3.1. 21 



The argii ment ation for this interpretation has been developed in lPaper 1 , Paper II 



IPaper Ilf and 'Paper V. The arguments will in the following sections be re- 
viewed and some points, which are less well documented in the papers will 
be substantiated. Experimental details given in the four papers will not be 
repeated here. 
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Figure 3.14: Reproducibility test in connection with iPaper Vl A 3D map was 
obtained 40 times of the same part of reciprocal space around a 400 reflection 
from a Cu sample deformed in tension to 1.8% strain. Each part is a projection 
onto the azimuthal plane, with qx horizontal and <jz vertical covering respec- 
tively 0.14A^^ and 0.006A^^. The datasets shown are separated by ~ 13min 
each. 



3.8.1 Observations from raw data 



The first high resolution data was prese nted in iPaper E and later a more de- 



tailed data set was presented in Paper III The observations presented are gen- 



eral for all Cu and Al samples investigated. 
The following data are from the experiment described in Paper ij . 



A copper sample was deformed to 2% strain, and a grain with a (100) direc- 
tion close to the tensile direction was identified using the "Reflection based" 
technique described in section 13.3.11 The grain was found to be smaller than 
the beam, both by spatial scanning and a calculation based on the full scattered 
intensity. The grain was hence fully illuminated during the measurement. A 
full map of the reciprocal space around one 400 reflection, consisting of 300 
a;-slices, was obtained. The rocking interval was chosen to be 0.007° corre- 
sponding to a size in reciprocal space of Aqz = 8 • 10~^A^^ hence close to the 
resolution in the other directions. An example of 9 a;-slices is shown in figure 
l3.2| on page [31] 

As the broadening is mainly in the azimuthal directions {qx and ^z)) it is feasi- 
ble to make a projection of the full 3D intensity distribution onto this plane (for 
details about this projection see section |3. 6. 2| |. The result of such a projection 
is a 2D map of the azimuthal spread, where the degree of freedom of strain 
has been integrated out. This projection is shown, for the full dataset, in figure 
13.151 With the contour levels used the map comprises isolated peaks (corre- 
sponding to the ones seen in the raw data), large islands and a rather smooth 
cloud of intensity, well above the detector background level. In the insert it is 
seen that the islands are composed of individual peaks, some of them partly 
overlapping. 
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Figure 3.15: Projection of the 3D reciprocal space intensity distribution onto 
the azimuthal plane ((Jx, (Jz)- The color scale of the contour lines represents 
integrated intensity in units of detector counts. The insert is a zoom on the high 
intensity part around {qxr^jz) = (—0.01A~^, —0.06A~^) with different contour 
lines and color scale. The circles indicate the position of 14 indiv idual p eaks 
selected for a statistical analysis as presented is section |3.8.2| From lPaper II . 



3.8.2 Single peaks 



The properties of the indivi dual pea ks were initially reported in iPaper 11 and 
analyzed in greater detail in iPaper II . 



14 individual peaks were selected for detailed analysis from the dataset pre- 
sented in figure I3ll5l 



The peaks were analyzed in the manner presented in section l3 .6 . 3l obtaining the 
following properties: Position in reciprocal space, full width at half maximum 
of the peak in the three directions, Aq = {Aqx, ^qij, ^qz) and integrated inten- 
sity of the peak. The integrated intensity was converted to a volume by the 
procedure described in section [3.6.51 Figure l3.16l show the width and volume 
data for the peaks. 
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(a) The fitted full width at half maximum, (b) Volume. The relative uncertainty is esti- 
Shown are the results for the three direc- mated to be ~ 15%, whereas the uncertainty 
tions in reciprocal space Aijj (x), Aqy (o), on the absolute scale of the volume is 50%. 
Aijz (□). The full horizontal line indicates The last two data points are from peaks in 
the average value of all peaks in all direc- the islands of enhanced intensity, 
tions {Aq) f» The error bars 

are calculated from the uncertainty in the 
level of the diffuse cloud. The fitted value 
is in all cases an upper limit for the width in 
the ^,y-direction. 



Figure 3.16: Properties of the 1 4 individ ual peak as shown on figure l3lT5l as a 
function of peak number. From lPaper 11 . 



Volume 

The mean volume from the 12 first, well-separated, peaks are {V) = 1.4p.m'', 
corresponding to an equivalent sphere diameter of (d) ~ 1.4(a.m. The width 
of the last two peaks is below the average width of the full set of peaks al- 
though the volume is larger. This indicates that the peaks do not arise from 
overlapping peaks from different regions. They do either come from larger 
dislocation-free regions, or they are from a part of the grain which is closer to 
the center of the vertical beam profile. 



Coinpared to TEM i nves tigations on similar polycrystalline systems (e.g. llEssmann , 
Il963t [Hansen et all 12001 1 and figure 11. 3t it can be seen that the length scale 
found above is comparable to what is reported for the dislocation-depleted re- 
gions in the dislocation structure (the subgrains). 



Quantitative cell size measurements exist on single Cu crystals iGottler , 1973 1. 



Estimating the equivalent single crystal resolved shear stress by correcting for 
the grain size strengtheni ng, it is p ossible to compare our polycrystalline data 
to the single crystal data |Paper lt \. It is found that the subgrain size in single 
crystals deformed to the resolved shear stress obtained here is ~ 2|^m. This 
is in good agreement with what we observe. Furthermore, it is kno wn that 



the d istribution of sizes of the dislocation-free regions are wide (e.g. llHuang , 
199811 ), consistent with the observation of larger scattering regions. 
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Dislocation density 

The mean width of the peaks in all three directions is (Aq) ~ 1.9 • lO^'^A^^, not 
far from the instrumental resolution as reported in section [331 

A large dislocation density would give rise to a substantial broadening of the 
peaks, this indicates that the dislocation density is low in the scattering regions; 
consistent with the interpretation as a diffraction signal from the subgrains. 



Spatial correlation 

The connection between the individual peaks and the spatial p osition in the 



grain was investigated as part of the data obtained for iPaper II . A number of 
(small) reciprocal space maps were obtained at a fixed tensile strain of 4.07%, 
while translating the sample with respect to the beam. The translation was by 
steps of 2|^m both in the x- and j/-direction. For this particular experiment the 
beam size was 14(a.m by 14|^m. 

In figure 13.171 projections onto the azimuthal plane are shown, as function of 
displacement from the position where the grain nominally is centered with 
respect to the beam. 

The important thing to notice is that each peak seems to be associated with one 
spatial position in the grain as the intensity in all cases seem to describe the 
intensity behavior of a bell-shaped beam profile. The beam is assumed to be 
bell shaped in both the horizontal and vertical directions due to the small size 
(see figure I3l4t . Information on the spatial po sition of the scattering entities . 
can i n principle, be derived from such scans [ PoulserJ, 20041: Lauridsen et al[ 

llool. 

That the peaks seem to be associated with only one spatial position in the grain 
strengthens the interpretation as a signal from subgrains. 



3.8.3 The diffuse cloud 



The diffuse cloud was in Paper ll tentatively interpreted as coming from the 



dislocation-rich walls in the dislocation structure. This interpretation was based 
on the fact that a large dislocation density will give rise to a broad strain and 
orientation distribution, hence to a signal in the reflection as the observed 
cloud. This interpretation has been substantiated firstly by an analysis of the 
mean strain in the material giving rise to the cloud [Paper 111] and secondly by 
measurements on a system having a dislocation structure similar to the walls 



I Paper V|] . 



3.8. DETAILED ARGUMENTS FOR THE INTERPRETATION 



63 




Figure 3.17: Azimuthal projections of 3D maps obtained at different positions 
with respect to the center of the grain. Each projection covers horizontally 
0.08A^^ and vert ically 0. 02A~^. Full information on the experimental condi- 
tions are given in lPaper 1 . 
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Figure 3.18: Position of the maximum intensity {cjy ^max), along lines in the qy- 
direction through the 3D intensity distribution, as a function of maximum in- 
tensity (/max)- Each average qy,max marked by a circle (Q) is taken over a small 
interval in intensity. The square (□) indicates the average qy,max position of 
all maxima with an intensity higher than 10000. Corresponding filled symbols 
indicate the number of points in each intensity bin. The horizontal full and 
dash-dotted lines indicate the maximum and average of the integrated radial 
peak profile, respectively. From .Paper IIL 



Statistical arialysis of interisity distributiori 



The statistical analy sis metho d presented in section 13.6.41 was applied to the 
data set presented in lPaper IIll The grain analyzed is from the same Cu sample 
and same conditions (2% tensile deformation) as the one investigated in the 
previous sections. 

Figure |3. 1 8 1 shows the results of the analysis. 

What can be observed is that a strong correlation exists between the maximum 
intensity and the radial position of the maximum. The dislocation-rich walls 
in a d islocation structure are, according to the composite model of Mughrabil 



1 1983ll , subjected to a forward stress (relative to the mean stress of the grain). 



and the dislocation-free subgrains to a back stress. If the model is true, and the 
interpretation of the cloud is true, a correlation as the one observed, would be 



the consequence llPaper IIlll 



This is not a direct proof of the correlation between the walls and the cloud, 
but as other results shown that subgrains are subjected to a backwards strain 
see section 14721 and fiPaper lEL it significantly strengthens the hypothesis. 
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Figure 3.19: Transmission electron micrographs on a grain in a 10% deformed 
AlMg sample. The grain is oriented 9° from the [100] orientation. Scale 
bars represents 200nm. A) The developed high density of dislocations. B) 
High ma gnification view, shows a imiform distribution of dislocations. From 
Paper Vl 



Non-cell-forming metals 

Measurements on a non-cell-forming AlMg alloy (aluminum + 4% magne- 
sium) were performed at d ifferent s trains (experimental procedures and fur- 



ther results are presented in lPaper V ). 



By TEM it was found that the dislocation structure, even after 10% tensile de- 
formation, was rather homogeneous. Figure |3 . 1 91 shows examples of electron 
micrographs. Figure l3.20l shows a comparison of raw a;-slices obtained on the 
AlMg sample and a Cu sample strained equivalently. No distinct peaks are 
seen in the data from the AlMg sample. Projections onto the azimuthal plane 
of the ob tained reciprocal space maps (as shown in figure 13721] and figure 4 in 
IPaper V ) likewise show no distinct peaks, but some smooth large-scale inten- 



sity variations exists. 

This indicates that a random tangle of dislocation gives a smooth intensity dis- 
tribution similar to the cloud observed from cell-forming materials. Further- 
more it shows that no clear peaks are observed when no subgrarns are present, 
substantiating the interpretation of the peaks and cloud observed for copper. 
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Figure 3.20: Raw a;-slices. Part A: AlMg, Part B: Cu. Cu data are from the 
dataset presented as Grain I in section l43l (a lso prese nted in lPaper ll) and the 
AlMg data are from the dataset presented in lPaper vl (see also table lOJ. Read 
from left to right the images show consecutive a;-slices. The images cover the 
same number of pixels but due to the difference in scattering and rocking angle 
the angular ranges covered is different, they cover 0.49° (AlMg) and 0.42° (Cu) 
in the azimuthal (horizontal) direction, 0.041° (AlMg) and 0.039° (Cu) in the 
radial (vertical) dire ction an d is integrated over (rocking angle) 0.003° (AlMg) 
0.0067° (Cu). From iPaperVl 
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Figure 3.21: Projection of a partial 3D reciprocal space map of the 400 reflection 
fr om an Al Mg sample deformed to 4.2% in tension. From the data presented 
in lPaper vi 
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3.9 Comparison to other techniques 

High Angular Resolution 3DXRD can be compared to other techniques on a 
number of parameters. The technique is fundamentally a 3D reciprocal space 
mapping technique, and should be compared to equivalent methods. On a 
general level the technique has to be compared to other techniques which are 
able to measure elastic strains and orientations of individual subgrains in a 
deformation structure. 

A advantage, that the technique inherits from the general 3DXRD principle, 
is the possibility for investigating individual grains deeply in the bulk of a 
polycrystalline sample. It may be debated to what extend this is important for 
the investigation of deformation structures. The use of polycrystalline samples 
has the advantage that grains of multiple orientations can be investigated in 
the same sample, and that the grain size limits the illuminated volume along 
the beam. This eliminates the need for complications such as the conical slit 
[Nielse n et al, 2000] or wire scanning [.Larson et ah, 2002]. 

The penetration power, due to the high energy, is clearly an advantage over 
traditional X-ray techniques. Copper Ka X-ray radiation has an attenuation 
length of w 22|am in copper whereas the here used X-ray radiation with an 
energy of 52keV has an attenuation length of 0.47mm. This ensures that true 
bulk informations are gathered. 



Reciprocal space mapping 



Traditional 3D recip rocal space maps are acquired point-by-point us ing an an- 
alyzer crystal (e.g. iFewster and Andrew , 119951: iMurphy et ali l200lll ). This al- 
lows for a better resolu tion in some recipro cal space directions, e.g. along the 
diffraction vector as in iMurphy et al\ [|200l|[ . The disadvantage of such tech- 
niques is that 3 degrees of freedom have to be scanned to obtain a full 3D re- 
ciprocal space map. 

The use of a 2D CCD detector clearly is an advantage with respect to acquisi- 
tion time, as only one degree of freedom has to be scanned. The disadvantage 
of using a CCD detector is that its point spread function will influence the re- 
ciprocal space resolution. The severeness of the effect is determined by the 
sample-to-detector distance (as shown in equation l3.9at , and can be controlled 
by changing this. The contribution to the reciprocal space resolution from the 
detector is for the present setup of the same order of magnitude as the contri- 
butions from the beam properties (as can be seen in table 13. li on pagel47t. 



Strain measurements 



The elastic strain (or equivalent the lattice spacing) in the individual subgrains 
have previously been measured by e.g. the convergent-beam electron diffrac- 
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tion (CBED) technique (e.g. llKassner et a/.Ll2002 !l), and ca n also be measured by 
the monochromatic version of the 3D crystal microscope iLevine et aZ.Ll2006ll . 

The instrumental resolution in the (jy-direction given by equation l3 . 14l (pagel46l 
corresponds to a width of 1.4-2.3 -10"^ in strain. However, the position of the 
maximum of a peak from an individual subgrain ca n be dete rmined better than 
this. We estimate the uncertainty to be ±1/2 pixel jPaper Iltl , corresponding to 
an uncertainty on strain difference of 0.7 • 10^*. 

The uncerta inty on the s t rain m easurements on the 3D crystal microscope as 
reported in llLevine a/. , 20061 is of sim ilar magnitude. The CBED measure- 



ments reported in | Kassner et ah , 2002|l have an uncertainty of ±1 • 10 ^ for 



aluminum and ±4 • 10 * for copper. 

This clearly shows that the X-ray based techniques have a better strain resolu- 
tion. Further more. High Angular Resolution 3DXRD has the advantage that 
problems with stress relaxations near the surface are eliminated by measuring 
in the bulk. 



Subgrain rotation 

The width of the instrumental resolution in the (j^-direction (as reported in 
equation 13.141 corresponds to an spread in lattice plane orientations of 0.006° 
- 0.008°. This is again an upper limit on the observable orientation differences. 
In the other azimuthal direction much better resolution can be obtained on 
expense of acquisition time by decreasing the rocking interval size. 

Orientation differences between different part s of a grain can be measu red by 
TEM techniques to an accuracy of 0.3° - 0.1° jLiul, Il995t IkuangL |2006|| . High 



Angular Resolution 3DXRD provides, at least, an order of magnitude better 
angular resolution. 



Summary 

These comparisons show that the developed technique, on the selected param- 
eters, has as good or better resolution as other existing techniques with respect 
to probing the properties of individual subgrains in a deformation structure. 

The main disadvantage of the technique is that no spatial information are ob- 
tained. Spatial information could be obtained by spatial scanning techniques, 
but it would, as always, be at the expense of data acquisition time. 

The main advantage, on the other hand, is the ability to perform relatively fast 
experiments on bulk structures. No other technique presently allows for inves- 
tigation of the d5mamics of individual bulk subgrains during deformation. 



Chapter 4 

Results and discussion 



4.1 Introduction 



High Angular Resolution 3DXRD was used for a number of experiments. The 
results provide new insight into the static properties of subgrains, and on the 
dynamics of the dislocation structures during deformation. Table \4A\ provides 
an overview of the experiments presented, and the sample materials used. 

The main results have been reported in six papers, as indicated below: 



Paper j presents the general technique and interpretation. The main scientific 



focus is on the d5rnamics of the subgrain structure during deformation 
(subgrain refinement), but a number of other issues are briefly touched 
upon. 



Paper 11^ gives a thorough description of the technique and the single peak 
analysis method. The static properties of the subgrains, namely the dis- 
tribution of elastic strain and dislocation density, are analyzed. 



Paper II j presents the statistical analysis method developed. The interpreta- 



tion of the cloud is substantiated, and the average strain in the walls is 
discussed. 



Paper IVl summarizes the main results of the first three papers. The focus is on 
communicating the results in a less technical and more materials science- 
focused framework. 



Paper Vl presents measurements on a non-cell-forming alloy (AlMg). This study 
is in this thesis mainly used as evidence for the interpretation of the ob- 
served cloud (as presented in section |3. 8.31 1. 
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Paper V j presents results from measurements during continuous deformation. 
The mea suremen ts complement the results on subgrain djmamics pre- 
sented in lPaperl , and provide independent results on the stability of the 



subgrain structure during stress relaxation and offloading. 

The results will be summarized in this chapter, collecting some of the frag- 
mented results which have been separated over more than one paper. A num- 
ber of results not yet included in any paper will also be discussed. Not all 
experimental details described in the papers will be repeated in this chapter. 



The two first sections (|4.2I and I4.3|l deal with t he static p roperties of the dis- 
location structure. The results from Paper I to Paper ivl are combined. Two 



main issues are dealt with: the distribution of elastic strain in the dislocation 
structure, and the dislocation density in the subgrains. 

Sections 14.41 and 14.51 deal with the dynamic properties of the dislocation struc- 
ture prob ed by in-sihi deformation experiments, the results from IPaper ll and 
Paper Vl l are combined. The first section (section 14. 4^ discusses the formation 



of subgrains during deformation of a fully recrystallized sample, and the possi- 
ble consequences of ending a deformation. The stability of the overall structure 
is furthermore analyzed with respect to stress relaxation and off-loading. The 
process of subgrain refinement is discussed in the second section (section |43J. 





Strain/Stress 


Maps 


Sample 


Discussed in 


Static (Grain I) 


Static 2% — 92MPa 


Full axial 400 


Cul 


Paper II 


Static (Grain II) 


Same sample and conditions as grain I 


Full axial 400 


Cul 


Paper III 


Static (Grain III) 


Static 2% — 82MPa 


Partial axial and side 400 


Cu II 


Continuous from 0% 


Strain rate 6 • lO^'^s^^ 


Small axial 400, Atjz w 9 • lO'^A"! 


Cu II 




Continuous pre-deformed 


Initial strain 1.82%. Two loadings at 1.1 • lO^^s"^. 
One loading at 3 • lO^^s^^. Final strain 3.137% 


Small axial 400, Aijz w 9 • lO^^A"! 


Cu II 


Paper VI 


Offloading 


From 111.9MPa to 49.7MPa 


Small axial 400, Aqz w OMA"^ 


Cu II 


Paper VI 


Stepwise deformation 


3% to 4.2%; steps w 0.04% 


Small axial 400, Aqz ~ 0.04A"l 


Cul 


Paper I 


Non-cell-forming 


Static at 1.8%, 4.2% and 10% 


Same 400 reflection at all strains (partial maps) 


AlMg 


Paper V 






plus additional grains at 4.2% (small maps) 





All grains investigated have the tensile axis close (within w 10°) to a (100) crystallographic direction. 
Partial maps cover a substantial part of the azimuthal spread in the (jz-direction. 

A single map has been obtained for the static investigations, whereas many were obtained for studies of the dynamics. 

(a) List of the main data-sets used in this study 



Material 



Grain size 



Processing 



Gauge dimensions 



Cul 



Cu II 



AlMg 



99.99% pure OFHC Cu 



As for Cu I 



tion 

95% Al, 3.9% Mg, 0.29% Fe, w 34^tm 
0.53% Mn, 0.14% Si 



36ixm (ignoring twin bound- Cold rolled to a reduction of 80%, resulting 8mm x 3mm x 300i_im 
aries) thickness 300 |^m, fully recrystallized by an- 

nealing for 60 minutes at 500°C 
Many grains of ~ 30|j.m, very Cold rolled to a reduction of 80%, resulting 8mm x 3mm x 300|j.m 
inhomogeneous size distribu- thickness 300 ixm, fully recrystallized by an- 
nealing for 120 minutes at 450°C 

Cold rolled to a reduction of 90%, followed 8mm x 3mm x 500 |^m 
by a reduction of 50%, resulting thickness 
500 pin, fully recrystallized by annealing for 
60 minutes at 575°C 



Cu I and II are the same base material, processed at different times. 

(b) List of the sample types used for the experiments. 

Table 4.1: Overview of datasets and samples. 
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4.2 Distribution of elastic strain 

The radial broadening of the reflections from a plastically deformed crystal 
(e.g. a single grain within a polycrystalline sample), shows that some distribu- 
tion of microscopic elastic strain exists in the crystal (see section [2.3. 11 1. How- 
ever no information is obtained on the spatial distribution of strains in the crys- 
tal. The microscopic strains may e.g. be described by the same distribution in 
any region of the crystal, or differences may exist between different parts of 
the dislocation structure. An example of the latter is the differences between 
subgrains and walls as suggested by the composite model (see page [19). 

High Angular Resolution 3DXRD allows for a direct determination of the mean 
elastic strain in the individual subgrains, using the single peak analysis method 
(section |3. 6.31 1. The analysis does also gives a limit on the width of the elas- 
tic strain distribution within each individual subgrain. Furthermore the mean 
elastic strain in the wall material can be analyzed by the statistical analysis 
method presented in section 15.6.41 

The mai n p art of the results presented in this section has been reported in 
Paper III and iPaper III . 



Results 

The two samples investigated were deformed to 2% plastic strain in tension, 
from a fully recrystallized state. All investigations were performed under load, 
with a fixed grip spacing (details on stress and strain are provided in table l4.lt . 
The grains selected have the tensile axis close to (within 10°) a (100) orienta- 
tion. Grains were selected using the methods presented in section [331 This 
grain orientation ensures that the morphology of the deformation structure is 
a simple cell structure (see section [l.l.2t . 

The choice of grain orientation furthermore allow s for comparison to t he re- 
sults by Ungar an d co-workers [Unea r et al ]ll984a|BlMushrab i et al. 1986, 
Ungar et^zTll99lll . They report asymmetrical line profiles (see section [2.3. H 



from [001]-orientated single crystals deformed in tension. Their lowest degree 
of deformation is to a true tensile stress of 64.2MPa at a tensile strain of 2.8% 



I Mughrabi et fl/.l,ll986|l . This is not far from the 2% tensile deformation reported 
here, where for grains I and II we find a tensile stress of 67MPa when cor rect- 
ing for the grain size hardening by use of the Hall-Petch relation I Paper ij l. We 



will hence in the following discussion compare to their results. 

The datasets on grain I and II each consist of a full map of one 400 reflection. 
The reflections are from the {100} planes which are close to being perpendic- 
ular to the tensile axis (the axial case in the notation introduced in the section 
on "asymmetric line broadening" on page [19). The grains are from the same 
sample, and were found by the reflection-based technique (section |3.3.1t ; the 
full crystallographic orientation is hence not known. 
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Grain I (axial) 


Grain II (axial) 


Grain III (axial) 


Grain III (side) 




'/i/,max 




6.9480 


6.9482 


6.9470 


6.9579 








6.9473 


6.9470 


6.9447 


6.9588 


Ai 


+ A2[10-3A- 




6.4 


5.5 


9.4 


8.5 


Ai 


- A2[io-3A- 




0.8 


1.4 


2.8 


-0.9 



Table 4.2: Parameter for the integrated radial peak profiles, based on the 3D 
reciprocal space maps obtained under load for the three grains investigated. 
Ai and A2 are the half -widths of the lower and upper part of the profile respec- 
tively. Ai -|- A2 is the full width at half maximum, and Aj — A2 a measure of the 
asymmetry. Measurements on grains in the undeformed state suggest that the 
absolute position of the peaks from grain III might be slightly (~ 1.4 • 10^'^A^^) 
shifted to higher ^jy-values. 



The dataset on grain III is from a different sample investigated at a later beam- 
time. The full orientation of the era was found by the GRAINDEX method 
(see section |2.5.1| |, this allows one to obtain reciprocal space maps for more than 
one reflection from the grain. Two 400 reflections were investigated, one being 
the axial case as above, and one a side case (having the scattering vector, close 
to, perpendicular to the tensile direction). The reciprocal space maps were not 
covering the full azimuthal spread of the reflection. 



Integrated radial peak profile 

The average elastic strain distribution in the illuminated part of the grains, is 
analyzed by projecting the obtained 3D reciprocal space map onto the radial 
direction (see section l3.6.2t to obtain the integrated radial peak profile. Figure 
14. II shows such a profile for grain I. 

The profiles have been characterized by fitting to a split-pseudo-Voigt function 
(see section |3. 6. 2t , and some statistical calculations. The main parameters de- 
rived from the radial profiles are the mean value, the position of the maximum, 
the peak width and the asymmetry. The parameters for the individual grains 
investigated are reported in table l4.2l 

It may be observed that each profile is asymmetric. In the axial cases the in- 
tensity falls off more slowly to the low-ijy-side and in the side case the slow 
decay is ori the hig h-qy-side. This is in agreement with what is reported by 
lUngar et all ll984a| on single copper crystals deformed in tension. 



It is furthermore observed that the width, asymmetry and absolute position 
differ somewhat between the grains. This must be due to the inhomogeneities 
of the stress-fields at the grain level, especially for grain I and II which are from 
the same sample. 



^The tensile axis is in the [14.5 2.3 1] direction for grain HI. 
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Figure 4.1: The integrated radial peak profile for grain I. Projection of the 3D re- 
ciprocal space intensity distribution onto the radial direction (<jy). The full line 
is a fit to a split-pseudo-Voigt function. Vertical full and dash-dotted lines indi- 
cate the positi on of the maximum and the average of the profile, respectively. 
Adapted from lPaper ij . 



Elastic strain in and between subgrains 

The mean strains of 1 4 individ ual subgrains from grain I have b een analyz ed 
(results are reported in Paper II (as (jy-positions of maximum), and Paper iV (as 
elastic strain)). The ijy-position of the maximum of the individual peaks, are 
found by the single peak analysis methods presented in section 13.6.31 Figure 
14.21 shows the maximum (jy-positions of the individual peaks. Also shown is 
the position of the mean value and the maximum of the integrated radial peak 
profile for this grain. 

Three important conclusions can be drawn from this data and from the peak 



width of the individual peaks (as reported in figure 3.16(a) on pagel6T]|: 



The individual subgrains have different mean elastic strains. 

The distribution has a width of 2.9 • 10^^, if defined by twice the standard 

deviation of the peak positions. 

The mean value of the peak positions are shifted by a significant amount 
compared to the mean value of the integrated intensity profile. 
The difference corresponds to a reduction of the elastic strain with respect 
to the applied stress (a backwards strain) in the subgrains by 0.9 • 10^^, 
on average. 

The width of the elastic strain distribution inside the individual sub- 
grains has an upper limitH of 2.4 • 10^^. 



^This upper limit is calculated from an upper limit of 1.7 ■ 10^^ A^^ on the mean strain broad- 
ening of the individual peaks. This limit is obtained by correcting the observed mean peak width 
of 1.9 ■ 10^^ A^^ (as reported in connection to figure [3~16(a)) with the minimum radial instrumen- 
tal broadening as given in equation l3.14bl and a contribution from size broadening, in a quadratic 
form. 
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peak nr 



Figure 4.2: The fitted (jy^mlx* "^-positions for 14 isolated peaks as a function of 
peak number (same peaks as reported on figure lSTSl and lXTSt . The uncertainty 
in the peak position is conservatively estimated to be ±2.5 • lO^^A^^. The av- 
erage (jy-position of the individual peaks is {^'?y^ml"'"^ = 6.9479 ± O.OOOSA^^, 

with a standard deviation on the positions of~ I.O IO^^A^^. The dashed- 
dotted horizontal line indicates the mean value of the integrated radial profile 
(as shown in figure 14.111 and the full horizontal line the maximum of the inte- 
grated radial profile. Adapted from lPaper 111 




10 15 
Peak number 



5 10 
Peak number 



(a) Axial case. ((Jy'mir") 



■ 6.9472A-1 (b) Side case. (ijy"mlx'"") = 6.9577 



Figure 4.3: The fitted <7y,max -positions for single peak s from grain III (some 
of the peaks were identified by M. Prinz llPrinj, i2006ll ). The dashed-dotted 
horizontal line indicates the mean value of the integrated radial profile and the 
full horizontal line the maximum of the integrated radial profile (as given in 
table 14. 2] l. The dotted horizontal line indicates the average (jy-position of the 
individual peaks. 
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An equivalent analysis on grain III shows very similar results for both the axial 
and the side case. Figure 14.31 shows the positions of the individual peaks in 
both the axial and the side cases. The distribution of mean strain of the sub- 
grains is wide, with a mean value very close to the maximum of the integrated 
radial peak profile (compare table 14^2) 1. The width of the individual peaks is 



furthermore similar to the values reported in figure 3.16(a) 



Strain in subgrain boundaries 

The technique does not allow for direct access of the strain in the subgrain 
boundaries. However, it was in section |3 . 8 .3 l argued that the observed cloud of 
enhanced intensity can be attributed to the walls. 

The statistical analysis method (section [3.6.4t will give access to the mean strain 
in the walls if the correlation between wall material and the cloud is accepted. 
The method has been a pplied to g rain II and the result is presented in figure 



13.181 on page l64l (and in IPaper lis ). This analysis suggests that the walls are 



in-fact subjected to an enhanced elastic strain with respect to the mean strain 
in the grain (a forward strain), as the mean <jy-positions at very low intensities 
indeed are below the mean value of the fully-integrated peak profile. 

It should be mentioned that such an analysis requires measurements with very 
good counting statistics. Equivalent analyses on the other datasets show the 
same trend, but not as pronounced a backwards strain as the one observed in 
figure l3l8l 



Discussion 

The composite model by 'Mughrabil |l983ll predicts that strain differences ex- 
ist between the subgrains and the walls. The asymmetry of peak profiles ob- 
served from plastically deformed crystals has been rationalized by use of this 
model, and a "de composition scheme" has been developed I Mue hrabi et al[ 
19861: lUngar efdi Il984a ,b] (see also section [2.3. It . The asymmetric peak is in 



this framework divided into two parts; one from the subgrains and one from 
the walls (see figure l272l on page [191. These parts have then been analyzed 
independently with respect to e.g. dislocation density. 

The model/ technique has been used in a number of studies, but has also been 
criticized. The existence of internal strains has e.g. been debated I Kassner et al[ 
200llStraub et fl/.l,ll996ll . 



The present technique gives the possibility for testing some of the assumptions 
and predictions from the composite model in a direct way. 

Our results directly show that the predicted reduction in elastic strain in the 
subgrains exists on average. Similarly our measurements indicate that the en- 
hancement of elastic strain in the walls exists. 
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The composite model predicts that the backwards strain in the subgrain, as 
observed in the axial case, is accompanied by a forward strain in the side case, 
due to cross contraction. This leads to a reversal of the asymmetry of the radial 
peak profile, as for the one observed for grain III in table 14.21 The magnitude 
of the forward strain in the side cas e is related to the backwards strain in the 
axial case by the Possion ratio [Mu ghrabi et all 1986 1. The ratio between the 



/ subgrain \ 
\^y,max , 


> / 


(^r ) - 1" 


side 


7 subgrain \ 
V 'jy,max / 


/ 




axial 



two strains can be calculated for grain III as (using values from figure |43] and 
table SJ 



= -0.4. (4.1) 



The res ult is close to the Poisson ratio which is 0.42 for a (100) oriented single 
crystal iLandolt-BomsteinFi The observed ratio is seen to be consistent with 
the predictions of the composite model. 

The backwards strain in the subgrains is furthermore found to be of the order 
of 10 — 25% of the average elastic strain of the grain under load (as evaluated 
from the mean position of the integrated radial peak profile with respect to 
an undeformed sample). This corresponds well with the result by Ungar et al. 
where a back stress which is 1 2% of the av erage flow stress is reported for the 
sample deformed to 26.2MPa |Ungar et fl/.LlT9"91|] . 

It should be mentioned that 'Levrne et al! flOoS] recently have reported com- 
parable results based on measurements on the 3D crystal microscope (as de- 
scribed in section |2. 5. 2t . Their data are taken in the offloaded state, and only 
for the axial case. 

Even as our results indicate that the composite model is conceptually correct, 
they also suggest that the results of the decomposition scheme for asymmet- 
ric peaks have to be re-interpreted. The result of the decomposition scheme 
applied to asymmetric integrated radial peak profiles is, as mentioned earlier, 
two well behaved peaks, one attributed to the subgrains and one to the walls 
(see figure IZ2l . The analysis of the two peaks has been done applying tradi- 
tional line profile theory. This corresponds to the implicit assumption that the 
major contribution to the broadening and shape of the peak attributed to the 
subgrains is the strain distribution within the subgrains and not among them. 

Our results show that the distribution of strains between the subgrains is wider 
than the internal strain distribution in the subgrains. The analysis of grain I 
shows that the ratio of the two widths is 0.8. As the estimate of the strain distri- 
bution within the subgrains is a conservative upper limit this ratio is probably 
even smaller. 



'The single crystal Poisson ration is calculated as |si2/sii|, with the elastic compliances 
Si2 = -6.3(TPa)-i and s„ = 15.0(TPa)-i. 
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To summarize, we suggest the following interpretation of the two peaks ob- 
tained by the decomposition method: 

• The position of the peak corresponding to the subgrains represents the 
mean of the distribution of strains between the subgrains. 

A backward strain exists in a average sense 

• The width of the peak corresponding to the subgrains, is highly influ- 
enced (if not dominated) by the distribution of strain between the sub- 
grains, and hence does not represent the distribution of strains within 
the individual subgrains. 

The results in this section show that performing 3D reciprocal space mapping 
at a high resolution can provide very important results regarding the internal 
elastic strain distribution in cell-forming metals. Furthermore the results pro- 
vide a test of some of the fundamental assumptions in the interpretation of 
more traditional X-ray analysis such as line profile analysis. 



4.3 Dislocation density in the subgrains 

An other interesting parameter is the dislocation density inside the subgrains, 
p. The density of dislocations influences the flow stress of the individual sub- 
grains, as T o< .yjo, and hence needs to be understood for making flow stress 
models which take into account the in-homogeneity of the dislocation struc- 
ture. 

The decomposition scheme of asymmetric radial peak profiles, as discussed in 
the previous section and on page[l9j has been used for evaluating the disloca- 
tio n density in the subgrains and walls. Su ch measurements are e.g. reported 
in llUngar et al[ 1984a ': Mughrabi et al. , 1986 '! where the ratio of dislocation den- 



sity between the walls and subgrains is found to be of the order of 10. The 
findings however contradict some elect ron microscopy-based measurements 
f Argon and Haasen , 19931: Haasen , 1993 1, where it is found that the dislocation 



density in the subgrains are "many orders of magnitude" smaller than in the 
walls. 

Our results suggest, that the radial peak profile obtained from the composite 
model-based decomposition scheme, is in fact highly influenced by the distri- 
bution of strains between the subgrains. The peaks for the individual subgrains 
are substantially narrower, and hence indicate a lower dislocation density. 

The ideal way to evaluate the dislocation density would be to utilize an ap- 
propriate theory or model which connects the dislocation density and disloca- 
tion configurations to the resulting 3D diffraction peak (or some projection of 
this). However, this is not possible within the time-frame of this study, as it 
would require a complete characterization of the instrumental broadening of 
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the High angular Resolution 3DXED setup, allowing for deconvolution of the 
peaks obtained. Two approaches, av oiding thi s problem, will be discussed in 
the following (results also reported in lPaper^ . 

A discussion will be given as to what can be learned on the basis of the upper 
limit of 1.7 10-3A-1 on the mean radial width of the physical peaks (as was 
reported in the last section (see the footnote on page [74| . The lack of full infor- 
mation on the peak shape will be treated by elimrnatrng all parameters but the 
peak width from a line profile theory using literature values. 

The results suggest that the dislocation density is very low and hence that "sta- 
tistical" approaches might not be appropriate. A second approach has there- 
fore been devised. This approach is based on direct calculations of the peak 
shape from a given dislocation configuration in a subgrarn. Presently this anal- 
ysis is restricted to one dislocation in the center of a subgrain, but in theory it 
should be possible to generalize the approach. 



Density of redundant dislocations 

A classical way to obtain information on dislocation density from integrated 
radial peak profiles is through a model for the micro-strain-distribution from 
some dislocation distributions. Such a model is then combined with the gen- 
eral theory for the Fourier coefficients associated with the radial peak profile 
(as discussed in section |2.3. It . This allows for a fit of the parameters of the 
dislocation distribution, based on the measured peak profile. 

The result s by Ungar and co-workers [Ung ar et Ml |l984a''F ; Muehra bi et at 



ir986; .Ungar et al 



near et al.l 1199111 are obtained by use of the classical theory by Wilkens 
I Wilkena, 1970allc |. The Wilkens' theory is therefore used in the following anal- 
ysis, as values reported by Ungar and co-workers then can be used for the 
calculations. 

Wilkens introduces the concept of a restricted random dislocation distribution. 
A region where the dislocations follow a restricted random distribution can 
be divided into a number of sub-regions of equal size, with the same number 
of, randomly distributed, dislocations in each sub-region. The dislocations are 
not fully randomly distributed in the entire region, hence the name "restricted 
random". It has been shown that such a dislocation distribution avoids the 
logarithmic divergence of the energy with crystal size, wh ich is associated with 
a totally random dislocation distribution I WilkensI, Il970b 1 . 



The distribution is characterized by two parameters, the dislocation density, p, 
and the parameter M. Where M is proportional to the average number of dislo- 
cations participating in the stress-screening configurationqj. A certain mixture 
of dislocation types (characterized by Burgers vector and direction) is assumed 



*M = Rey/p with Rf the outer cut-of f radius, which is almost equal to the size of the sub- 
regions of the distribution IWilkensl , ll970ar . 
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when the Wilkens' theory is applied to the peak profile analysis. It is further 
assumed that the dislocation distribution, for each type of dislocation, consists 
of an equal amount of positive and negative Burgers vectors, hence the term 
"redundant dislocations". 

Our peak profiles do, as mentioned, not contain enough information to fit both 
parameters in a traditional Wilkens' analysis. We have therefore chosen to fix 
the M value. The Wilkens' theory then predicts the width of the integrated 
radial peak profile as function of the dislocation density. 

We have chosen to use the value of M = 1.60 given in |Ungar et al. , 1984a|| for 



a deformation which is comparable to ours. We furthermore assume the same 
mixture of dislocation types, allowing for use of the parameters describing the 
dislocation t5rpes given in lUngar et flUll984all . 



The limit of 1.7 • 10~^A~^ on the width of the physical peak profile, is for apeak 
which is not integrated in the azimuthal directions. Wilkens' theory predicts 
the shape of the azimuthally integrated radial peak profile. This possibly in- 
troduces an additional broadening of the peak, depending on the details of the 
3D peak shap^ Correction for this gives an upper limit on the azimuthally 
integrated radial peak profile of 2.6 • lO^'^A^^. 

Based on this upper limit on the width of the integrated radial peak the disloca- 
tion density can be calculated. The parameters by Ungar and co-workers and 
the theoretical connection between the normalized M factor and normalized 
half width given in figure 1 in [Wilkens, 1970a] are used in the calculation. The 
resulting limit on the dislocation density is 12 • 10^^m~^ s ubstantially lower 
that the value of 22 • lO^^m^^ reported by Ung ar et"a/1 | 1984a ]. It should be em- 



phasized that our result is a very conservative upper limit on the dislocation 
density, and that the real dislocation density possibly is much lower. 



Density of excess dislocations 

Excess dislocations (that are dislocations not compensated by a dislocation of 
opposite Burgers vecto r nearby) will h ave a tendency to broaden the reflection 
in the azimuthal plane IWilkens , 1984 1. 



Two arguments which show that as little as one excess dislocation in a sub- 
grain will lead to a visible "splitting" of the reflection, will be provided. Both 
arguments are based on considerations of one edge dislocation positioned in 
the center of a subgrain. 

In the first model the dislocation is treated as a part of a symmetrical tilt bound- 
ary consisting of edge dislocations with a separation of the subgrain size (D). 
Figure l4!4l' a) illustrates this situation. 



^The width of the integrated and non-integrated profiles differ at most by a factor of 1.55, 
which is in the case of a 3D peak profile following a squared Lorentz function. 
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a) b) 




Figure 4.4: Illustration of the two models of one dislocation in the center of a 
subgrain. a) The dislocation is considered as part of a symmetrical tilt bound- 
ary, b) The dislocation is considered as being placed in the center of a cylinder 
with free boundaries. 



The angular difference between the two sides of such a boundary, j6, is, as men- 
tioned in section ll.l.li given as /5 = b/D where b is the length of the Burgers 
vector (2.55A for (110)1/2 dislocations in the {111} plane for Cu). Setting 
D = 1.1 i^m (the mean 2D length scale of the subgrarns found in IPaper Il|) re- 
sults in /5 = 0.013°. If the two sides of the boundary can be considered as 
scattering, partly, incoherently this will directly lead to a splitting of the reflec- 
tion in the azunuthal direction of (?azimuthal ~ 

1.6-10-3A-1. 

The second model is based on the observation that subgrains with a volume 
of w l|^m'^ comprise ~ 3000^ Cu atoms. This limited size makes it feasible 
to make direct calculations of the diffracted intensity from a 2D section of a 
subgrain. The atomic positions in such a section can be obtained from the 
elastic displacement field from the dislocation configuration. 

The calculation is done on a cylindric fee crystal with a dislocation in the center 
and free boundaries. Figure fOf b) illustrates this configuration. The cylinder is 
infinite in the [100] -direction with a diameter of D in the perpendicular plane 
(comprising directions [010] and [001]). The crystal is orientated such that the 
crystallographic coordinate system and the experimental coordinate system co- 
incide. The edge dislocation in the center has a Burgers vector of [01 1] 1 /2. The 
direction of the dislocation line is the [100] -direction; hence it is not a mobile 
dislocation. With this geometry the position of the atoms can be calculated 
from the elastic approximation of the displacement field fo r one dislocation in a 
cylindrical crystal with free boundaries, as determined by iLeibfried and Luckd 
|l94 9]. The elastic and crystallographic parameters were chosen to be those of 

The reflections associated with such a system are essentially two-dimensional 
due to the translation symmetry of the crystal. The intensity distribution can 
therefore be calculated from the atomic positions in a single (100) plane; the 

^Poisson's ratio for polycrystaUine copper a = 0.34 and lattice spacing a 3.6149A^^ were 
used. 
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Figure 4.5: Atomistic simulation of the intensity distribution within the 040 
reflection for a quasi 2D fee crystal with one 1/2[011] edge dislocation in the 
center. In the calculation elastic and crystallographic properties of Cu were 
used and the boundaries taken as free. The cylindrical domain had a diameter 
of 3000 inter-atomic distances (w l.l|j.m) corresponding to an area of 0.9i.im^, 
giving a dislocation density of 1 • lO^^m^^. The gray-scale bar is in dimension- 
less units. 



calculation is 2D both in real and reciprocal space. The diffracted intensity 
distribution is obtained by standard kinematic diffraction theory as described 
in section lZTl using equation |2.3| and l2.4l (on page[T2t with the atomic scattering 
factor set to one. 

Figure 143] shows the resulting intensity distribution within the 040 reflection. 
It may be observed that the distribution splits into three parts each separated 
by w 1.3 • 10^3^-1 with the major splitting in the azimuthal t/z-direction. 

These two models of a single dislocation in a subgrain, show that High Angular 
Resolution 3DXRD might be sensitive to as little as one excess dislocation in 
a subgrain, as it will lead to a visible split of the associated peak. It should 
be noted that the orientation of the dislocation with respect to the scattering 
vector, determines the visibility of the dislocation in the diffracted signal. 



We have generally not found such split peaks in the static investigations llPaper Illl . 
indicating that the existence of single excess dislocations is rare. 

A stepwise loading experiment has been performed (see section |431 "stepwise 
deformation" in table l4.ll and Pa per I). The peaks from a number of subgrains 
were followed in-situ during stepwise straining. During this experiment we 
did observe several peaks which split into two at some strain and then recom- 
bine at a later strain step. Figure |4!6] shows such a sequence. The interpretation 
is that a single (or a few) dislocation(s) got trapped in the subgrain investi- 
gated, and that it /they at a later deformation stage got released. 
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Figure 4.6: Excerpts form larger azimuthal projections of reciprocal space maps 
taken at increasing strain. These are shown from left to right as a function of ex- 
ternal strain (with strain increments relative to the first map shown). The peak 
followed split into two parts (as indicated by the arrow) and then recombines. 
Colo r scale a nd contour lines are identical to those on figure 14.141 Adapted 
from Paper j . 



4.4 Formation and stability of subgrains 

Investigations of subgrains during continuous deformation are possible using 
High Angular Resolution 3DXRD. Such experiments make it possible to ad- 
dress questions about the formation and stability of the dislocation structures. 

Experimentally, the specifications of strain rate, wished time /stain resolution, 
acquisition time for one detector image (a;-slice), and the size of the 3D maps 
are highly rnterluiked. Typically the 3D maps are highly truncated in the qz- 
direction in order to obtain a reasonable time/ strain resolution. 

Two experimental scenarios have been used. In the simplest case; one focuses 
the beam on a particular grain, but does not try to keep the beam at a precise 
position within the grain. As function of time the same reflection from a grain 
can then be mapped, but not necessarily the same part of the reflection. In 
the more challenging case one tries to keep the beam fixed with respect to a 
position within the grain and perhaps to compensate for macroscopic grain 
rotations leading to movement of the reflection. 

In the following, two examples of such continuous deformation experiments 
will be presented. They both relate to studies of formation and stability of 
dislocation structures. 

A fundamental limitation of TEM is (as mentioned in section |1.1.3t that the 
technique is ex-situ (unless thin films are investigated). The consequences of 
offloading and sampl e preparation for TEM investigatioris hav e been consid- 
ered in the past (e.g. (Mughrabi', 1971 : \bung and Sherrilll-ll^ ). However, it 



is impossible to show by TEM if the dislocation structures observed after the 
deformation has ended are the same as the ones that existed during the defor- 
mation. 
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Fundamentally two very different scenarios for the consequences of terminat- 
ing the deformation can be imagined. The structural formation may be a part 
of the plastic deformation process, the dislocations generating the structures 
dynamically while the sample deforms, and hence nothing dramatic happens 
at the end of deformation. Alternatively the subgrain structure might be the 
result of a relaxation process at the moment where the deformation stops, and 
therefore not be representative of the dislocation distribution which existed 
during the deformation. 

It should be possible to distinguish between these two scenarios by continuous 
deformation experiments with the present technique. 

The first experiment addresses the initial formation of dislocation structures as 
a fully recrystallized sample is continuously deformed, but does not address 
the consequence of ending the deformation. 

The second experiment addresses the possible changes in a dislocation struc- 
ture at the moment where a continuous deformation is stopped. The influence 
on the overall dislocation structure by stress relaxation and offloading is finally 
discussed. 

4.4.1 Formation of dislocation structures 

By following a reflection from a grain while continuously deforming the sam- 
ple from the fully recrystallized state, it should be possible to observe the mo- 
ment where / if the reflection starts to break up, indicating the formation of a 
dislocation structure. 

This experiment has been tried a number of times; the fundamental experi- 
mental problem being that a grain has to be followed while it moves due to 
the straining. This is especially hard in the case of an undeformed sample, as 
the stress rig will deform (mainly elastically, but also some movements in gears 
and the like) in the initial part of the deformation. Furthermore, it might be dif- 
ficult to follow the reflection as the direction of the grain rotation is unknown. 
Two data sets which illustrate the general impression from these experiments 
are presented below. 

The first experiment was done with very few a;-slices per map, as a proof of 
concept (data not included in table 14.11 as the experimental conditions were 
somewhat different). The sample was tensile deformed from 0% to 3% with a 
strain rate of 2.5 • lO^^s^^. It was observed that the breakup did happen while 
straining, and that a clear breakup into sharp peaks superimposed on a clou d 
existed at all strains equal to or above 0.4% (this result is reported in lPaper I ). 

Additional data were taken at a later stage, in connection with the experi- 
ments reported in Pa per VL A grain with a size smaller than the beam size 
was centered in the beam and the corresponding 400 reflection centered on the 
detector. The sample was strained continuously from 0% with a strain rate 
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Figure 4.7: Stress-strain curve for the continuous deformation from 0% (see also 
table l4T| . The individual markers correspond to one a;-slice, and the colors 
indicate the strain ranges of the 3D maps. The arrows indicate the four maps 
illustrated in figure 14^81 The reason for the gaps (corresponding to missing 
datasets) is that the peak had to be re-centered due to grain rotations, and while 
this was done no 3D maps could be obtained. 



of 6 • 10^ s~ , while obtaining reciprocal space maps consisting of 15 a;-slices 
each. The (/z-rage covered was ~ 9 • lO^^A^^, and the center of the range was 
adjusted during the experiment to compensate for grain rotation. Figure 14.71 
shows the stress-strain curve for the experiment, indicating the individual o;- 
slices as markers, and the 3D maps with colors. It is seen that each map covers 
a small strain interval due to the continuous deformation. 

The projections of the 3D maps onto the azimuthal plane (using the pixel-based 
projection, as described in section |3.6.2t are shown in figure 14.81 for four suc- 
cessive maps. No breakups are seen before these maps, and the reflection was 
imfortimately lost after obtaining these maps, most likely due to grain rotation. 

In spite of the experimental problems, the two data sets give equivalent results. 
The breakup into subgrarns clearly happens during the deformation, and with- 
out any significant broadening of the peak before the point where the structure 
is seen. The breakups happen at very low degrees of plastic deformations, in 
the last case already at 0.12% deformation. 



4.4.2 Stability of dislocation structures 

A more controlled experiment can be performed by continuous deformation of 
a pre-strained sample. 

A grain with the tensile axis close to a (100) direction was identified, and the 
sample was pre-strained stepwise to 1.82% tensile deformation. A calibration 
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Figure 4.8: Projection onto the azimuthal plane of 3D reciprocal space maps 
obtained during continuous tensile deformation from 0%. The projections are 
shown as a function of stress and strain. The images span w 25 • lO^'^A^^ hor- 
izontally (corresponding to —qx) and ~ 9 • lO^'^A^^ vertically (corresponding 
to qz)- Note that the color scales are different at the different strain levels. 

curve for strain increment and vertical sample movement as a function of ten- 
sion motor movement was found by following the grain with t he X-ray b eam 
during this initial straining. Further information in table l4.1l and |Paper VI . 

The sample was thereafter strained continuously at different strain rates, and 
with different following stress relaxation times. Figure |4l9l shows the tension 
motor movement and the stress-strain curve for the full loading cycle, com- 
prising two slow loadings with strain rate 1.1 •10~^s~^, one fast loading with 
strain rate estimated to be 3 • 10~'^s~^, and the offloading of the sample. 

Small 3D reciprocal space maps, comprising 15 6t;-slices each, were obtained 
during the slow strarnrng and the holding periods. Figure 14.101 shows a se- 
quence of projections of such maps obtained during the second slow load- 
ing and holding. Larger maps were obtained before and after the offloading, 
projections of these maps are shown in figure 14.111 The 3D maps were also 
projected onto the radial direction and the integrated radial peak profiles ana- 
lyzed, the resulting peak positions and widths are shown in figure l4T2l 

What may be observed is that the subgrain structures develop and change as 
long as the deformation is in progress; this is seen directly from the azimuthal 
projections (figure l4.10t . The increase in the width of the integrated radial peak 
profile during the deformation, as seen on figure l4T2l indicates that the inter- 
nal strain distribution broadens as the sample deforms. When the deformation 
is stopped the overall structure seems to freeze, and only minor changes are 
observed in the structures as seen from the azimuthal projections. A clear re- 
laxation is however observed in the width and position of the radial profiles. 
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Figure 4.9: Loading cycle for the continuous deformation of a pre-strained sam- 
ple. Top) Displacement of tension motor. Bottom) Stress-strain curve. Verti- 
cal lines on the two figures identify the "events" of the loading cycle. From 
iPaper Vll 



The relaxation observed in the peak position is no surprise as this corresponds 
to the observed macroscopic decrease in stress. 

The visual impression of the minor changes in the azimuthal projection is that 
the individual peaks become sharper, corresponding to some cleanup pro- 
cesses in the structure, which is consistent with the decrease in width of the 
average internal strain distribution. 

This substantiates the conclusion from the last section, that is, the dislocation 
structural formation and development is part of the plastic deformation pro- 
cess, and the structures observed after stopping the deformation are the same 
as the ones existing while straining. The observed stress relaxation seems (on 
the timescale of the experiment « Ih) to be connected to a clean up process in 
the walls and subgrains, and not to major rearrangements of the structure. 

To further investigate when the subgrain structure starts to evolve, one sepa- 
rated peak was selected from the data obtained during the second slow load- 
ing. This case is special, as the sample had been left for stress relaxation for 
41 hours before the measurement, the stress relaxation was therefore substan- 
tial. The peak clearly rotates (shifts) with respect to the rest of the structures 
observed in the maps. Figure l4?T3l shows the macroscopic stress and strain, to- 
gether with the position of the mean (jy-position of the integrated radial peak 
profile and finally the i^i- -position of the selected peak. All data have been nor- 
malized to coincide at the start and end of the deformation. The time where 
the macroscopic stress reaches the previous maximum flow stress is further- 
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Figure 4.10: Projections onto the azimuthal plane of the last six 3D reciprocal 
space maps obtained during the second slow straining, and every second of 
the maps obtained after the straining was stopped. The time given is with 
respect to the time where the straining was stopped. The projections each cover 
0.12A~^, and 0.009A^^ in the and qz directions respectively. Adapted from 
IPaper vk 
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Figure 4.11: Projections onto the azimuthal plane of the 3D recipr ocal spa ce 
maps obtained before and after the offloading of the sample. From .Paper VI . 
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Figure 4.12: Position and width of the azimuthal integrated peak profiles as a 
function of time. The time is in each part measured with respect to the start 
of the deformation. 1. slow loading (open square), 1. holding (filled square), 
2. slow loading (open diamond), 2. holding (filled diamond), fast loading 
(triangles pointing up), offloading (circles). 
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Figure 4.13: Analysis of the correlation between the macroscopic stress (red 
line), macroscopic strain (blue line), mean elastic strain as found from the po- 
sition of the integrated peak profile ( x ), and subgrain peak position in the qx- 
direction (□). Data are taken during the second loading, after a long (41 hours) 
stress relaxation period. The vertical line indicates the point where the macro- 
scopic stress reaches the previous maximum flow stress. All four curves have 
been normalized to give best agreement at the start and end of the deformation. 
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more indicated. What can be observed is that the position of the mean value 
of the radial profile follows the macroscopic stress as expected from elastic ex- 
pansion. From the macroscopic stress curves it is furthermore seen that the 
material starts to yield before reaching the previous maximum flow stress, this 
must be due to the pronounced relaxation process which took place during the 
long holding period. The (j^-position of the peak is overall linearly dependent 
on the macroscopic strain, but a significant deviation from linearity is seen 
around the point where the macroscopic stress reaches the previous maximum 
flow stress. 

The deviation from the linear dependence between the (j^-position of the indi- 
vidual peak and the macroscopic strain starts a certain amount of strain after 
the material has started to deform plastically, but before the previous maxi- 
mum flow stress has been reached. This can be understood in the framework 
of the composite model, as it predicts a two stage yield process. Initially only 
the dislocation free subgrains will deform plastically and then later the harder 
walls also deform plastically. This suggests that the evolution of the subgrain 
structure only happens when the entire dislocation structure is deforming plas- 
tically. 

Before and after the offloading of the sample larger reciprocal space maps were 
obtained. Projections of these are shown in figure l4?TT] What can be seen is that 
the offloading of the sample does not give rise to changes in the overall struc- 
ture (some rotation is ho wever observed) . Th e result is in line with a recent X- 
ray based experiment by lSchafler et all j2005|] , who report in-situ synchrotron- 
based line broadening studies of [100] orientated single crystals. At selected 
strains Schafler and co-workers interrupted the deformation, unloaded the 
crystal, reloaded and continued the deformation lUngan, 1200611 . Peak profiles 



were obtained at all stages, and no differences existed between the average dis- 
location density during loading, under load, and in the offloaded state tUngail, 
2006ll . 



Beside the investigation of the overall structural changes due to offloading 
seven individual peaks were analyzed using the single peak analysis meth- 
ods presented in section 15.6.31 The peaks clearly move in the azimuthal direc- 
tions during the offloading. The analysis shows that the movement is mainly 
due to elastic rotation of the full structure (average shifts of — 0.0072A^^ and 
— 0.0024A^^ were found in the qx- and <jz-direction respectively). However 
a significant spread was observed in the individual rotations, the standard 
derivations being 13% and 20% of the average rotation in the qx and qz di- 
rections respectively. 

The (/y-positions of the individual peaks showed that a backwards strain (as 
discussed in section |4!2|| also exists in this case. What may be observed is that 
the magnitude of the mean backwards strain increased b y 14% upon the of - 
floading. This is consistent with the results presented in liBorbely et ali l2000ll , 
which are based on traditional line profile analysis on Cu single crystals de- 
formed under constant load. More interestingly the width of the distribution 
of strains between the subgrains decreased by the significant amount of 31%. 

This shows that the general structure is independent of offloading. However, 
offloading does change the internal strain distribution and the exact orientation 
differences between the subgrains. 
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4.5 Subgrain refinement 

A prominent feature of the d5mamics of a dislocation structure is the cell refine- 
ment process, that is the fact that the average size of the subgrains goes down 
as the plastic deformation increases. 

Traditionally this refinement is attributed to the trapping of dislocations inside 
subgrains, which leads to the buildup of new boundar ies with cor respond- 
ing orientation differences between the two sides (e.g. iHughes et ali Il997ll ). 
The signature of such a process would be that the peak from the initial sub- 
grain started to broaden, as dislocations get trapped, and then divided into 
two peaks rotating away from each other. 

It is possible to follow such a process by mapping the same region in reciprocal 
space as a function of strain. The results and experimental details are also 
reported in lPaperl . 



A grain, having the tensile axis close to a (100) orientation, was identified. 
The sample was pre-strained to 3% elongation and then from 3% to 4.2% in 
steps of ~ 0.04% while following the grain. The beam used was reduced to 
14|^m X 14|^mby slits. The illuminated part of the grain was kept constant with 
respect to the center of the grain by spatial scanning of the grain at (almost) all 
strain steps. This scanning procedure has the consequence that possible beam 
drift problems (as described in section 15. 2. It are corrected for. 

Figure l4.14l shows a projection onto the azimuthal plane for one of the 3D re- 
ciprocal space maps obtained. By keeping the stain increments small and the 
maps relatively large (in the (jz-direction), it is ensured that the possible rota- 
tion of the individual subgrains can be followed in detail. By making movies, 
consisting of projections onto the azimuthal plane as a function of strain, it is 
verified that the individual peaks can be followed (o ne such movie is published 
as "Supporting Online Material" in connection with lPaper j) . 

An increase in strain from 3% to 4.2% corresponds to an increase in the number 
of subgrains of w 66% assumi ng an inverse squa re root law of the boundary 
spacing as a function of strain I Hughes et a/.Lll997ll . By this is should be highly 



likely to observe the signature of subgrain refinement, as more than half the 
subgrains have to be "refined". 

Approximately 20 peaks were followed for some strain intervals, and in no 
cases were a clear signature of a subgrain breakup observed. However, we did 
observe an unexpected intermittent dynamics, peaks appearing out of regions 
of enhanced intensity, exist for some strain steps before they disappear again. 
Figure |4. 151 shows an example of such an event. Our interpretation of this is, 
that the subgrains "materialize" out of regions of high dislocation density, stay 
for some time and then "dissolve" again. 

The picture of intermittent dynamics constitutes an alternative framework for 
understanding the d5mamics of dislocation structures. It may e.g. explain how 
a preferred orientation of dislocation boundaries are kept during cold working. 
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Figure 4.14: Projection on to the azimuthal plane of a 3D reciprocal space map 
obtai ned at 3 .49% tensile strain from the stepwise deformation experiment. 
From Paper j . 
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Figure 4.15: Excerpts form larger azimuthal projections of reciprocal space 
maps taken at stepwise increasing strain. These are shown from left to right 
as a function of external strain. Strain increments are relative to the first map 
shown. Example of a peak appearing, existing for some strain interval (note 
the ■ ■ ■ indicting that multiple strain steps where the peak was visible, have 
been left out), and then disappearing. Color scale and contour lines identical 
as in figure 14141 Adapted from Pa per L 



Chapter 5 

Conclusions and outlook 



The main aim of this study has been in-situ investigations of bulk deformation 
structures in simple polycrystalline fee metals at low degrees of deformation. 
Both the dynamics of such structures during deformation and the static prop- 
erties under load have been investigated. 

The outcome has been two fold; a novel technique has been established and a 
niraiber of fimdamental scientific questions have been addressed. 

The technique 

The novel s5mchrotron radiation-based technique "High Angular Resolution 
3DXRD", which has been developed, extends the 3DXRD microscopy tech- 
nique to high resolution 3D reciprocal space mapping. The technique has 
a number of advantages over traditional 3D reciprocal space mapping tech- 
niques: 

• Individual broadened reflections from deeply embedded single grains in 
a polycrystalline sample can be mapped, due to the high energy (52keV). 

• The mapping is relatively fast due to the high flux and the fact that only 
one degree of freedom has to be swept to obtain a 3D map. The typical 
acquisition time for one a;-slice is 10s. This enables in-situ experiments 
during deformation, as a stress rig has been incorporated in the setup. 

• The voxels in reciprocal space, mapped onto one data point, are close 
to being cube-shaped, giving an equivalent resolution in all directions. 
The dimensions of the voxels are w 1 • 10~^A~^ in the radial and in one 
azimuthal direction, if defined by the instrumental broadening. The di- 
mension in the perpendicular azimuthal direction is defined by the rock- 
ing angle interval. 
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The resolution of the technique corresponds to a strain resolution of 0.7 • 10~ 
or better. This is as good as or better than other existing techniques for mea- 
suring the strain in individual subgrains. 

The orientation resolution is « 0.007°, which is an order of magnitude better 
than TEM-based methods. 



Interpretation of 3D reciprocal space maps from cell-forming metals 

The technique has been used to investigate polycrystalline copper at low de- 
grees of plastic tensile deformation (< 5% strain). 

It was foimd that the broadened reflections, from such a cell-forming metal, 
comprise a structure consisting of bright sharp peaks (separable in 3D in the 
reciprocal space) superimposed on a cloud of enhanced intensity. Based on the 
size of the scattering entities, the width of the peaks and spatial separation it 
is concluded that the peaks arise from individual subgrains (dislocation-free 
regions) deeply embedded in the grain investigated. The cloud is inferred to 
arise from the dislocation-filled walls in the dislocation structure. 

The technique developed provides a unique direct probe of the properties of 
individual subgrains in bulk grains, something which has been unavailable 
previously. A number of scientific issues regarding the dislocation structure 
have been investigated using this technique. 

Static results 

Based on static investigations of grains, with the tensile axis close to a (100) 
orientation, in deformed samples under load it is fotmd that: 

• The subgrains on average are subjected to a reduction of the elastic strain, 
compared to the mean strain of the full grain, when lattice planes perpen- 
dicular to the tensile axis are investigated (the axial case). The size of this 
strain reduction is on the order of 10-25% of the average elastic strain in 
the grain. An increase in the elastic strain is by contrast observed from 
the lattice planes parallel to the tensile axis, and the ratio of these two 
"internal" strains is close to the Possion ratio. 

• Evidence imply that the walls in the dislocation structure are subjected 
to an increase in elastic strain. 

• The distribution of elastic strains between the subgrains is broader than 
the distribution of elastic strains within the individual subgrains. 

• The redimdant dislocation density in the subgrains is very low (< 12 • lO^^m"^). 
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• Calculations show that as few as one excess dislocation in a subgrain 
might give a visible splitting of the associated diffraction peak. Such tran- 
sient splittings have been observed diiring stepwise loading. 

The observed elastic strain distribution suggest that the composite model by 
Mughrabi correctly describes the elastic strain distribution on average, but that 
asjnrraietric peak profile analysis has to be reinterpreted. 



Dynamic results 

A number of in-situ experiments have been performed, investigating the prop- 
erties of the dislocation structure during and after deformation. 

The development of the dislocation structure was followed during continuous 
deformation from the undeformed state. The subgrains clearly form during 
the deformation, and seemed to start forming as soon as plastic deformation is 
detectable (in one case investigated at around 0.12% elongation). Additionally 
a pre-deformed sample was continuously deformed, monitoring the subgrain 
dynamics diiring the deformation and immediately after termination of the de- 
formation. The evolution of the dislocation structure clearly happened duiing 
the deformation, and no changes were observed at the termination point. 

The consequences of w Ih stress relaxation after continuous deformation and 
of imloading of the sample were also analyzed. It was found that the over- 
all dislocation structure only depends on the maximum applied stress, as no 
overall changes were observed in the subgrain structure during the relaxation 
(at least on this time scale), or after the offload. However, a minor relaxation 
in the width of the internal strain distribution is observed from the width of 
the integrated radial peak profile. On the level of the individual subgrains, it 
was observed that the width of the stain distribution between the subgrains 
decreases by a substantial amoimt (31% in the case investigated) and that the 
level of backwards strain increases slightly (« 14% in the case investigated) 
during the offloading. 

Based on stepwise loading of a pre-strained sample it is found that the dislo- 
cation structure shows "intermittent" d5mamics. The subgrains seem to arise 
out of regions with an enhanced dislocation density, exist for some time with 
proceeding deformation, and then eventually disappear. In none of the investi- 
gated cases was a simple subgrain breakup observed. The subgrain refinement 
process hens does seem to be due to this intermittent djmamics. 

The results of the dynamical experiments points to a dislocation structure which 
is very d5mamic. The structure seems to adjust itself to fit with the present en- 
vironment by creating and dissolving subgrains. Furthermore, it seems that 
major changes happen only when the entire dislocation structure imdergoes 
plastic deformation. 

In conclusion it is my belief that measurements of this kind provide very valu- 
able information on a scale and under conditions previously unattainable. Hope- 
fully, the results can inspire a new generation of work hardening and pattern 
formation models. 
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CHAPTER 5. CONCL USIONS AND OUTLOOK 



Outlook 



The results presented in this thesis have shown that High Angular Resolution 
3DXRD is a very promising technique for investigating subgrains and their 
d5mamics. 

One of the most surprising results is the observation of intermittent dynam- 
ics. However subgrain refinement experiments at higher strains are needed to 
show if these findings are general for plastic deformation o r only apply at the 
early stages of the pattern formation (see e.g. comment by KubinI 1 2006t] ^. 



Furthermore a number of additional scientific cases, which are in line with the 
present study, are candidates for being studied by High Angular Resolution 
3DXRD: 



Strain path change. As the deformation structure depends on the mode of de- 
formation, it will need to adapt itself if the strain path is changed. It 
would be very interesting, and possible, to investigate if the structure 
gradually changes from one to another by e.g. rotation of subgrains, or if 
the change is mediated through a process of "dissolving-and-generation" 
as with the case of intermittent dynamics observed in unidirectional de- 
formation. 

Grain orientation dependence of internal strain. The grains investigated in this 
thesis all have an orientation which results in a cell-like morphology of 
the dislocation structure. The possible grain orientation dependency of 
the internal strain distribution can be studied by investigating grains of 
different orientations, and mapping multiple reflections from each grain 
(as was done with low resolution in [ Pantleon et ah, 2004]). Such experi- 
ments can e.g. show if the observed forward and backwards strains also 
exists in grain with the more general cell-block structure. 

Investigations of highly deformed materials. The investigation of highly de- 
formed metals is a region of material science which is attracting a lot of 
attention presently. We have performed proof of concept experiments on 
thin films and one a;-slice from such a measurement is shown in figure 
15.11 It is seen that it is possible to observe individual spots from the sub- 
grains from such a highly deformed sample. Hence it should be possible 
to investigate e.g. internal strain differences in such materials. 

To be able to fully perform such experiments a few experimental developments 
are needed, to improve the statistics and stability of the measurements: 

Automated peak finder for 3D reciprocal space maps. The individual peaks are 
presently located manually from azimuthal projections of 3D maps. This 
is rather time consuming, and it is not feasible to investigate all peaks in 
a large map. An automated procedure would give better statistics, and a 
more unbiased peak selection. 
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Figure 5.1: A single a;-slice (covering 0.01°) taken on the 220 ring from a thin 
film of highly deformed aluminum. Two clear peaks are seen on a rather 
smooth background. The degree of deformation is unfortunately unknown 
but very high (von Mises strain > 5). This is from the very first experiment, 
therefore the 3DXRD peak shape analysis setup was used (see section |2.5.3t . 



Full 3D fitting of peaks. The individual peaks are presently analyzed by fit- 
ting line profiles through the peak, the full 3D information is therefore 
not used in the analysis. An automated 3D fitting routine should be able 
to give better information on the integrated intensity, and the true shape 
of the peak. Furthermore 3D multi peak fitting may allow for better sep- 
aration of partially overlapping peaks, and hence give access to more 
peaks. 

Beam position monitor and feedback systems for improved stability. The ma- 
jor experimental problem is presently the instability of the beam position. 
This has the consequence that the relative intensity between peaks might 
change over time, due to beam drift, which can only be partly corrected 
for by tedious spatial scanning of a grain with respect to the beam. A 
feedback system should be able to keep the position of the beam con- 
stant, and thereby the intensity profile on the sample constant. 

Even higher resolution and/or better characterization of the instrumental resolution. 

The instrumental broadening contributes significantly to the width and 
shape of the individual peaks from the subgrains. If precise measure- 
ments of internal strain distribution or dislocation density in the individ- 
ual subgrains are wished for, it is important to be able to fully separate 
the contribution from the instrument and the real physical peak profile. 
This could be obtained in two ways; firstly if the instrumental broaden- 
ing could be characterized better, a deconvolution could be performed, 
secondly the instrumental broadening could be reduced so much that its 
contribution is negligible. The major problem with characterizing the 
instrument is that it is hard to find a probe which can be used for the 
measurement. If small grains in a powder are used, size broadening will 
influence the measurement, and if a perfect single-crystal is used, full 
d5rnamical scattering theory is needed. Reduction of the instrumental 
broadening is not easy either as all contributions are of the same order of 
magnitude, hence major parts of the setup would have to be improved to 
give a significant improvement. However, a detector with a smaller point 
spread would improve on this (a reduction by a factor of 2 in the point 
spread would give a reduction of the theoretical experimental broaden- 
ing by 20%). It should be mentioned that the Darwin width of a Cu 400 
reflection is w 3 • lO^^A^^ only a factor 3-4 away from the present reso- 
lution. 
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One important observation is the diversity of the 3D shapes of the peaks from 
the individual subgrains, which varies from flat disks through spheres to nee- 
dle shapes. This must be because of differences in dislocation configurations 
and boundary conditions. In section |43] results using traditional line profile 
analysis and simple atomistic calculations were presented. It should be possi- 
ble to obtain more information from the 3D shape of the peaks. 

A possibility might be to make full 3D simulations of subgrains with embed- 
ded dislocations and given boundary conditions. Such simulations could give 
an insight into how dislocation configuration and density influence the shape 
of the 3D diffraction peak. Ultimately it might be possible to make a fitting 
routine which matches the calculated peak to the measured 3D peak using the 
dislocation configuration as the "fitting parameter". 

Finally I emphasize that. High Angular Resolution 3DXRD as such is not lim- 
ited to investigations of dislocation structures and their dynamics under defor- 
mation. Examples of other applications are: 



General elastic strain distributions. An equivalent investigation on strain dis- 
tributions to the one presented in section 14.21 can be performed on any 
sample which consists of a number of incoherently scattering entities (e.g. 
crystallites in a powder). By 3D reciprocal space mapping it will be possi- 
ble to investigate if each entity has an internal strain distribution or if the 
strain is mainly distributed between them (a previous example of such 
results exists in jpewster and Andrew , Il999ll on the grain scale in a poly- 



crystalline sample). Such investigations could be of major importance for 
interpretation of e.g. powder diffraction signals. 

Peak separation. A general condition for 3DXRD experiments is that a separa- 
tion of individual diffraction peaks is needed. High Angular Resolution 
3DXRD might be a possible route for separating reflections which overlap 
when observed by the (relatively) low angular resolution of the normal 
3DXRD microscope. 

Annealing experiments. Annealing phenomena in the deformation structure 
can be investigated by including a furnace in the setup. This could e.g. 
be for investigations of the recovery process of highly deformed samples. 
Such experiments would benefit from the ability to separate individual 
peaks, and wo uld allow for investig ating not only the volume and ori- 
entation (as in fGundl ach et al[ '2004] where traditional 3DXRD is used 
for such experiments) but also allow for monitoring the internal strain 
distribution as a function of time. 
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